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- Kleven Year Variation of Cosmic-Ray Disturbance and Its 
Relation to Solar and Geomagnetic Activities 


By Yukio MIYAZAKI and Masami WADA 


Scientific Research Institute 


Abstract 


The daily mean intensity of cosmic rays and geomagnetic 

horizontal component are measured from the 27-day running average 
values, and this measure of disturbance is named “ variability.” 
The frequency distribution of variability for both cosmic rays and 
horizontal component deviate from the normal distribution toward 
negative, and this character is nearly the same amount for two 
cases throughout one sunspot cycle from 1936 to 1945. 

The maximum value of cosmic-ray variability delays from the 
sunspot maximum. Since the circumstance is nearly similar to that 
of horizontal component, it may be said that a large number of the 
cosmic-ray disturbance observed at the sunspot maximum time is due 
to the phenomena related to the magnetic disturbances. 

Though there are several evidences which indicate the close 
relation between cosmic-ray and geomagnetic variability, apparent 
difference is seen at the period of sunspot minimum. It suggests 
that there is some kind of magnetic disturbance which is independent 
of cosmic-ray variation. 

1. Introduction 

It has been well known that the intensity of cosmic rays changes at the time 
of magnetic storm [1], and that the intensity itself shows the 27-day recurrence 
during several intervals as was established by Monk and Compton [2]. While 
magnetic activity also shows the 27-day variation, the problem still remains what 
relationship does exist between these two kinds of variations [3]. Recently, Meyer 
and Simpson [4] have reported that the amplitude of the cosmic-ray 27-day intensity 
variation from 1936 to 1953 is closely related to the eleven year cycle in general 
solar activity. Similar problems have been studied by the present authors from 
somewhat different point of view. The results are presented in this paper and 
possible interpretations are considered. 

The authors aimed not only at the cosmic-ray 27-day variation, but also at 
any disturbance having period of several days, hence, we carried on the difference 
of the daily average from the 27-day running average intensity. This measure 
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may be called the variability of cosmic-ray intensity (abbrev. V-/). Since magnetic 
horizontal component is thought to be closely connected to the change in intensity 
of cosmic rays, the same procedure mentioned above was applied on it (V-H). 
For the variability of cosmic-ray intensity, the Forbush ion chamber data for the 
period 1936 through 1946 at Huancayo (0.6° south geomagnetic latitude) [5] were 
analysed, and for the magnetic variability, Kakioka Magnetic Observatory data in 
Japan were used. Kakioka (26.0° north geomag. lat.) is situated at comparatively 
magnetically less disturbed latitude than the polar region and hence the day-to-day 
variation of the horizontal component is thought to indicate the field having the 
close relation to cosmic-ray variation. 

Making use of the variability one can disregard the variation in periods 
longer than 27 days. However, because of the 27-day running average, apparent 
increases of the intensity are seen before and after the actual severe decreases. 
Such distortion restricts the method of analysis. Therefore, we applied some 

_ statistical treatments on the variabilities. 

The character figure of dark Ha- flocculi which was used in the analysis of 
M-region by Waldmeier [6], as well as relative sunspot-numbers for the year 1936 
through 1946 were used as the measure of the solar activity. 

2. The frequency distributions of variabilities 
The V-Jand V-H histograms are shown in Fig. 1. We can easily notice that 
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Fig. 1. Histograms of the cosmic-ray (full line) and geomgnetic (broken 
line) variabilities summarized for years from 1936 to 1945. 
Nnmber of samples is 3400. 


they display the same distribution. Fig. 2 is the quantitative representation of these 
distributions for each year. The lower part of Fig. 2 is the degree of asymmetry 
which is the cubic root of g,-statistics [7] derived from the third moment around 
the mean value. If the distribution is normal, this value comes to zero, and owing 
to the finite number of samples, 95 percent of observed points would fall in the 
hatched region in the figure. The degree of sharpness, which is the quadruple root 
of the absolute value of g.-statistics derived from the fourth moment around the 
mean value, is plotted in the upper part of Fig. 2. The hatched region has the 
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Eleven Year Variation of Cusmic-Ray Disturbance 


same meaning as above. Fig. 2 indicates 
that the distribution of V-J and V-H 
in each year is not normal, but is largely 
deviated towards negative value. No 
regard was taken of the differene among 
years, because the scattering of value 
may be of little meaning. 

The frequency distribution of 
cosmic-ray variability at the interna- 
tional quiet or disturbed days are 
shown in Fig. 3. It is seen that the 
large deviation towards negative is 
mainly due to the magnetic storm effect. 

The standard deviation of each 
year is shown in Fig. 4 together with 
the annual mean of sunspot-numbers, 
in which the scale values are normalized 
by the mean of the whole period. From 
this, one may find that V-J runs almost 
parallel with sunspot-numbers, and this 
result is just the same as the amplitude 
of semi-annual value of the 27-day varia- 
tion by Meyer and Simpson [4]. Never- 
theless, the maximum position of V-J 
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Fig. 2. The degree of asymmetry (lower part) and 
the degreee of sharpness (upper) obtained 
froin the histogram of each year. Hatched 
region is the error of sampling from the 
population of normal distribution with 95 
percent fiducial probability. 


appeares a year later than that of sunspot-numbers, and the circumstance is nearly 


similar to V-H curve. The range 
of the variation of sunspot- 
numbers in the eleven year 
cycle is the largest of three. 
In the case of V-I, we 
subtracted the natural fluctua- 
tion. Without this correction, 
the results are almost the same. 
3. The 27-day mean deviations 


SS of the variabilities 
¥-I PERCENT The mean deviation which 
Fig. 3. Histograms of the cosmic-ray variability is 0.8 of the standard deviation 
corresponding to the international quiet in th 
(broken line) and disturbed (full line) days. im the case of normal distribu- 


Number of samples is 550, 


tion is obtained by averaging 


the absolute value of variability for the Bartels’ 27-day solar rotation period. They 


are plotted in Fig. 5 together with the 27- 


day average of sunspot-numbers and 


the monthly value of dark Ho-flocculi character figures from Mt. Wilson Observatory 
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reported in “ Terrestrial Magnetism and Atmospheric Electricity.” 

There are several peaks in V-J and V-H curves. V-I peaks gradually 
decrease towards the sunspot minimum, while those of V-H do not show any change 
with sunspot-numbers but the peak widths gradually become wider. In spite of 
these difference in the form of peaks, V-I and V-H peaks almost coincide with 
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each other, from 1936 
through 1945, 16 peaks coincide 
out of 35 during 120 solar rota- 
tions. Now, 120 rotations were 
divided into three parts accord- 
ing to the grade of sunspot- 
numbers: (a) period of sunspot 
maximum, (b) period of decrease, 
and (c) period of minimum. 
The results are shown in Table 
I in which the product of J and 
H means the number of coinci- 
dent peaks between V-J and 
V-H, while that of J and H 
means the number of V-J peaks 
not accompanied by V-H peaks 


i.e, 


and vise versa. 

The probability of having 
no relation between V-J and 
V-H is 8% in the case of (a) 
and 0.4% in (b), hence, we may 
conclude that there exist the 
relation between V-J and V-H. 
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Relation between the peaks (in Fig. 
5) of 27-day mean deviation of cosmic- 
ray and geomagnetic variabilities. 
White circles are those in the period 
from 1936 to 1945, and black circles 


are in 1946, 


Fig. 4. 
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The standard deviations of the cosmic-ray 
(black circle) and geomagnetic (white circle) 
variabilities and sunspot-numbers in each 
year for 1936 through 1945. The ordinate 
scale values are normalized by mean of 
whole period. 


Table I 


ee Gee (c) ) 


| | H sum H | H sum H H sum 


a. Bf Sah AR a Re 10 12 


ae 
I Dl bs a 6 | 23/29 7 21 | 28 
% sum 12 28 40 | 14 | 26 40 | 9 oa | 40° 


On the anne ‘the i in the cese of 
(c) is not significant. 

Fig. 6 is the relations between V-J and 
V-H of coupled peaks. There seems little 
evidence that V-J and V-H are in simple 
relation. This scattered quantitative relation 
is the cause of the fact that the significant 
level of (a) in Table I is not so much higher 
than that of (b), but lower. There are fre- 
quent occurrence of disturbances in the period 
of (a), and they overlap with each other in 
the 27-day period which make the time of 
peak rather obscure. 
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Though the analyses presented here are unpreferable methods for the study 
of 27-day recurrence phenomenon, we can derive the following informations from 
them. In Fig. 7 which is drawn by superposing the peaks in Fig. 5 for V-F and 
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Fig. 7. Superposed average curves of 27-day mean deviations of cosmic-ray 
and geomagnetic variabilities. 35 peaks in Fig. 5 (all in 1936-1945) 
are taken as original days. 


V-H, respectively, a peak appears only at the original point, and no decay curve is 
seen. It may be said that the amplitude of recurrence phenomenon is small when 
it appears again, so that they cannot be discriminated from the background 
phenomena which may be the after effect of disturbances of such kind. 

However, broad peaks of V-H curve at the time of sunspot decrease in Fig. 
5 are seemed to be the appearance of 27-day recurrence of magnetic activity called 
as M-region [8]. Meanwhile, V-J peaks are not clear at that time and moreover, 
they drop to the order of the natural fluctuation at the time of sunspot minimum. 
From these facts, we may conclude that the fractions which still remain in V-H 
at sunspot minimum are not effective to the cosmic-ray disturbances. 

Meyer and Simpson [4] delated Forbush type variation from their study on 
the amplitude of 27-day variation. The periods which were delated from are shown 
in Fig. 5 by arrows. There might still exist some large peaks in their study, and 
their result may be somewhat ambiguous in this respect. We do not commit 
ourselves to the problem of 27-day variation, because the 27-day variation of cosmic- 
ray intensity and the Forbush type magnetic storm effect have not yet been clearly 
distinguished from each other. 

One may notice that there is difference in background lines between V-I 
and V-H in Fig. 5. Namely, eleven year variation is more remarkable in V-J than 
V-H. However, we shall study about it in next paper. 

Character figure of dark Ha-flocculi shows clear eleven year variation, but it 
does not reach to zero at the sunspot minimum period. This is because we choose 
character figure instead of the area. If one takes area, such deviation above zero 
may be vanished [9]. Restricting the curves in Fig. 5, we can not expect for 
different effect from different measures of solar activity. 
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4. Discussion and summary ; 

The variations of the disturbances of cosmic rays and geomagnetism in the 
sunspot cycle have been examined above. Here we intend to summarize the 
known facts, and to discuss by what degree these facts might be established or 
not. 

Cosmic-ray intensity variation has been observed at the time of magnetic 
storm. This effect is also recognized by the peak coincidence of V-Jand V-H, and 
at the same time, by the fact that the frequency distributions of the variabilities 
of them show almost the same character. In spite of this, any simple quantitative 
relation have not been found from coincident peaks, and it is the already known 
fact for the case of individual storms [1]. 

The difference in V-J and V-H at the time of sunspot minimum suggests that 
at least two types of magnetic disturbances exist. Recently, Kitamura [10] reported 
that remarkable decrease of cosmic-ray intensity was observed at the time of SC 
storm, but there was little effect at non-SC storm, even if the disturbance index 
was large. Meanwhile, Newton and Milson [11] have studied the relation between 
sunspot-numbers and the occurrence frequency of SC and non-SC storms, and 
concluded that SC storm has closer relation to sunspot-numbers than non-SC storm, 
and the latter has stonger tendency of 27-day recurrence than the former. These 
results do not seem to contradict with the fact that the occurrence of cosmic-ray 
disturbance is almost parallel to the sunspot-numbers. 

On the other hand, the strong evidence of 27-day recurrence tendency of 
cosmic-ray intensity which was established by Monk and Compton [2] is not agreeable 
with the phenomena mentioned above. We wonder whether this phenomenon is 
due to SC storm or non-SC storm effect, and whether it is independent of magnetic 
disturbance or not. 

Finally, we may summarize the resluts of the analyses as follows: 

(1) The frequency distributions of deviated values in daily mean cosmic-ray 
and geomagnetic intensities are of the same character as seen in Figs. 1 and 2. 

(2) The eleven year cycle variations are seen in the disturbances of cosmic 
rays and magnetic force, and the phases nearly coincide with each other. 

(3) However, the times of their maxima are not identical with the sunspot 
maximum, but shift several years. In view of this fact, disturbance of cosmic-ray 
intensity relates closer to geomagnetic activity than sunspot-numbers. 

(4) Magnetic disturbance does not become small at the time of sunspot 
minimum, but cosmic-ray disturbance becomes very small even to the amount of 
natural fluctuation. This fact suggests that the cosmic-ray variation relates more 
simply to solar activity than geomagnetic disturbance does. 

(5) Further, it may be considered that the magnetic disturbance is caused 
by other component which is quite independent of cosmic-ray disturbance at the 
time of sunspot minimum. 


In conclusion, the authors would like to express their thanks to Drs. S.E. 
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Forbush and I. Lange for the publication of cosmic-ray data on which the present 
research was performed. 
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Magnetic Interaction between Ferromagnetic Materials 
Contained in Rocks* 


By Seiya UYEDA 
Geophysical Institute, Tokyo University 


Abstract 


The detailed mechanism of the self-reversal of thermo-rema- 
nent magnetism of rocks was studied. The ferromagnetic mineral 
grains responsible for this phenomenon were determined as A’B- 
grains, which consist of two constituents, A’ and B. The magnetic 
and crystallographic examinations of Haruna rocks showed that the 
constituent with higher Curie point, in the two-constituent model, 
should be estimated to be a solid solution between ilmenite and 
hematite, and not the Ti-poor titanomagnetite as was estimated to 
be in the previous studies. In other words, the magnetic interaction 
causing the self-reversal of thermo-remanent magnetism is between 
two ilmenite-hematite solid solutions with different Curie points. 
Electron-microscopic observations showed that the constituent with 
higher Curie point is intergrown into the constituent with lower 
Curie point as fine lamellae. The observed configurations of these 
two constituents indicated that the phenomenon of reverse thermo- 
remanent magnetism could be explained theoretically, provided that 
the grains as fine as the observed lamellae have very strong 
thermo-remanence. 
Similar examinations were conducted on ferromagnetic miner- 
als of several other rocks with the same mineral assemblage as the 
Haruna ferromagnetic minerals. The results of these examinations 
indicated that there exist, in natural ferromagnetic minerals, various 
degrees of magnetic interaction which are governed by their state 
of co-existence. 
§ 1. Introduction 
In the previous papers [1], [2], [3], written jointly by Nagata, Akimoto and 
this author, the discovery of the self-reversal phenomenon of thermo-remanent 
magnetism of igneous rocks as well as the tentative results of the study on its 
mechanism were reported. In outline summary those results were as follows: the 
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ferromagnetic minerals contained in a dacitic pumice of Mt. Haruna and in a 
dacitic pitchstone of Mt. Asio mainly consist of two phases. One is a titanomagnetite 
phase, called A, and the other is a phase of ferromagnetic solid solutions of 
ilmenite and hematite, called B. Their Curie points are respectively about 500°C 
and 230°C. In the process of a weak field-cooling from a temperature above the 
Curie point of B, the thermo-remanent magnetism of B is produced in the direction 
just opposite to that of the external magnetic field, due to a magnetic interaction 
with the magnetization of A, which has been magnetized in the direction of the 
external field while the temperature of the specimen was still above the Curie 
point of B. 

Although the composition of the ferromagnetic minerals is the same for 
both the Haruna and the Asio rocks, there is an important difference in their 
characteristics: i.e. while the total thermo-remanent magnetism of the Haruna 
rock is reversed, the Asio rock shows the reversed magnetization only in some 
‘cases of partial thermo-remanent magnetism. This difference, though essential 
from the geophysical point of view, was interpreted to have been caused merely 
by the fact that the amount of reversely oriented magnetization of B in the Asio 
rock was not large enough to overcome the normal magnetization of A at the 
room temperature in the total thermo-remanent magnetism. 

To clarify the actual mechanism of the above mentioned interaction, some 
detailed study was made on the various magnetic behaviours of individual single 
grains of ferromagnetic minerals in these rocks. According to that study, the 
reverse thermo-remanent magnetism takes place within almost all the individual 
grains containing B, while the grains containing only A have the normal charac- 
teristics of thermo-remanence alone. This fact suggests that the mechanism of 
reverse thermo-remanent magnetism in the present cases is somewhat different 
from Néel’s original model [4]. 

Since the discovery of this peculiar phenomenon, many scientists in palaeo- 
magnetism examined the possibility of the occurence of similar phenomenon on a 
large number of the naturally reversely magnetized rocks, in an attempt to explain 
their reverse natural remanent magnetism [5], [6], [7], [8], 9], [10]. But in most of 
these tests, the results were reported to be negative. These negative results have 
promoted the wider acceptance of the hypothesis that the earth’s magnetic field was 
reversed at the times of formation of these rocks. However, in a few cases among 
them, it was reported that the magnetic interaction between two different con- 
stituents should not be overlooked. For instance, Balsley and Buddington [11] have 
shown that the reverse natural remanent magnetization of the rocks in Adirondack, 
U.S.A. would be explained by the two-constituent mechanism; Graham [12] also 
mentioned the possibility of some cases where the magnetic interaction may play 
an important role in magnetization of rocks; Saito [13] discovered recently another 
example of partial self-reversal of thermo-remanent magnetism in an iron sand at 
the Sokota-mines, Japan, which is quite similar to the Asio case. 
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As can be expected, magnetic interaction seems to exist always when magnc- 
tically different phases are intimately intergrown, whatever the phases of the 
concerned ferromagnetic minerals may be. In the present paper, special attention is 
given to the study of the general nature of the magnetic interaction acting in rocks 
containing the ferromagnetic solid solutions of ilmenite and hematite, together with 
the titanomagnetite phase, placing special emphasis on the Haruna-type self-reversal. 

Abbreviations of technical terms used in this paper are as follows: 


gi bl <d% Lara ustees thermo-remanent magnetism, 

RTRM ......reverse thermo-remanent magnetism, 
H.,..........externally applied magnetic field, 

Je Menace Sete demagnetizing field, 

TERME: . : titanomagnetite, 

Eten. os. Ss ilmenite-hematite. 


§ 2. Nature of Grains Responsible for Production of Reverse Thermo-remanent 
Magnetism 

As was stated in the Introduction, the self-reversal of TRM of the Haruna 
type was found in a dacitic pumice of Mt. Haruna, a pitchstone of Mt. Asio and 
an iron sand of the Sokota-mines. Since the experiments [3], [13] revealed that 
their self-reversal mechanism must be quite similar, the Haruna rock will mainly 
be discussed in this paper. 

During the study of RTRM characteristics of single grains [3-2], the peculiar 
fact was noticed that intense RTRM can take place even in the aB-grain, where 
the presence of Ti-poor Ti-Mt phase was hardly detectable by thermo-magnetic 
analysis. As long as the two-constituent mechanism theory is accepted, this fact 
may force one to choose one or the other of the following two possibilities : (i) the Ti- 
poor Ti-Mt phase playing an essential role in RTRM can be assumed to be so scarce 
that it is less than the limit of detectability by the thermo-magnetic study, say 1%, 
or (ii) the effective A-constituent in the two-constituent model is not Ti-poor Ti-Mt 
but some other phase that has a J,(T)-curve hardly distinguishable from that of B 
when they are mixed. At the stage of the previous study, it was more natural to 
assume the first proposition. However, a number of later experimental results 
have made the second proposition more plausible. 

Fig. l,a is a photomicrograph of the polished surface of a grain showing a 
distinct stepwise character in the J,(7)-curve. In this figure, the fact that coarse 
eres of A and B are attached to each other can be observed. A simple calcula- 
tion can show that such a state of coarse co-existence has little relevance in 
ei Sees the appearance of RTRM. On the other hand, on the surface of the 
grain of which the J,(T)-curve showed only the Il-Ht type curve, no Ti-Mt phase 
was observed (Fig. 1,b). Both of these grains have the RTRM characteristics. 

’ It a . concluded, from these observations, that the presence of the A 

phase attached to the B phase in the AB-grains i iti 

the appearance of RTRM and that bretere = pear: spas, te 
uld be excluded from 
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Fig. 1,a. Photomicrograph cf an AB-grain, Fig. 1,b. Photomicrograph of an aB-grain, 


showing coexistence of A- and B- showing;no existence of A-phase in 
phases in coarse masses. Reflected visible size. Reflected _ light. 
light. 


the study of its mechanism. 

In order to exclude the trivial A phase, a chemical method was used together 
with thermo-magnetic separation. Since HCl effectively dissolves Ti-Mt but not 
ferromagnetic II-Ht, thermo-magnetic separation at about 250°C and boiling in HCl 
for 5 hours of the whole specimen of the Haruna ferromagnetic minerals were 
alternately repeated for the elimination of the Ti-Mt phase. The non-existence of 
Ti-Mt phase in the grains thus extracted was ascertained by the fact that they 
gave no tinge of Fe***-ions to HCl. From 206gr. of the original ferromagnetic grains, 
7.26gr. of such grains was obtained by the above method. Therefore, the ratio of 
the latter to the former 0.035 may be regarded as giving the approximate value of 
the concentration of grains responsible for the RTRM in the original ferromagnetic 
minerals in Haruna rock. In Fig. 2, the full curve shows the mode of the develop- 


ment of the total TRM of the grains thus obtained, while the broken line indicates 


that of the original sample. In comparing these two curves, the fact that the above 


separated grains are actually responsible for the RTRM will be understood. 

It will also be seen that, along the field-cooling curve of the separated grains 
in the figure, there is a distinct rise of the value of normal magnetization at about 
210°C below which the curve goes steeply downward. This rise, that was not 
distinct in the curve of the original sample, may be considered to be due to the 
very magnetization that exerts H, upon B when the sample is cooled through the 
Curie point of B. The curve will, then, suggest that the Curie point of the A 
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phase responsible for the appearance en ee aie 

of RTRM is about 320°C and not as 

high as that of the Ti-poor Ti-Mt, as 

was presumed to be the case in the 
previous study. Thus a new 50 
constituent was found which more 
adequately explains the RTRM than 

the Ti-poor Ti-Mt constituent, in the 
two-constituent model. In order to 0 
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distinguish these constituents, the 
new one will be called A’. Therefore 
the grains obtained by the above stated 
separation will be called A’B-grains. 
As for the phase of this A’-constituent, ~5° 
there are two possibilities: a Ti-rich 
Ti-Mt phase, or an_ II-Ht solid 
solution, like the JB-constituent, 
but with a higher Curie point than  igg 
the latter. An X-ray analysis of the 


A’B-grain offered some _ indication ee ; IT = 
Sai ig: 2. velopment of total TRM: J(T)="(T)x Hex 
that the second possibility shoul r 

P lipid +Jrun,,(T), where H..=2.0 Oe. Full 


regarded as the more probable, as line: A’B-grains. Broken line: original 


will be stated in § 4. ferromagnetic mineral grains. 
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the estimated thermo-magnetic Curie point is shown also in larger scale 
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ak constituent. H,, by broken line, H,,=22,00¢, 
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The results of a careful measurement of the temperature change of spon- 
taneous magnetization and of susceptibility, namely the J,(7) and the «(T) of the 
A’B-grains are shown in Figs. 3,a and 3,b, where the curve for the Jatter. is 
plotted also in larger scale near the Curie point. The kink at about 260°C, though 
not very distinct, may have been caused by the included A’. If the above argument 
that the relevant A-constituent is a ferromagnetic material with Curie point at about 
320°C, or A’, is correct, the fact that RTRM takes place most remarkably during the 
field-cooling in the limited temperature interval of 300°C and 250°C in the original 
specimen can be explained very plausibly, since this temperature interval is 
supposed to include the domain-fixing temperature of A’. In the previous paper 
[1-4], where the Curie point of the relevant A-constituent was supposed to be about 
500°C, this characteristic was explained by the assumption that the reversible part 


of the magnetization of the Ti-poor Ti-Mt would be primarily responsible for H, 
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Fig. 4,a. Modes of development of total TRM and Fig. 4,b. Field dependence of Jr(T))=*(T») 
its decay by heating in non-magnetic space, xH,, and of total TRM 
for various intensities of H,,. (A’B-rgains) J%  (T). (A'B-grains). 
Full lines: modes of development: J(7) 350, Her 
=K(T)XHeetJq y,,(T)- Broken lines: 
modes of decay of fixed magnetization: 


KT)=S 7° 4, (7). 
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exerting on B. 

The A’-constituent, being supposed to be intergrown in the host of B, may be 
considered to have been generated by an exsolution process. Consequently, the 
chemical composition of both A’ and B may depene on the local density of the 
exsolution in the host. This means that the Curie points of A’ and B may vary 
over a certain temperature interval. The Curie points of A’ and B estimated 
from the above expriments would, then, be the mean values. 

The modes of development of total TRM and its demagnetization by heating 
in a non-magnetic space for various field intensities are shown in Fig. 4,a, and the 
«(T))xH,, and the i - (T,) curves deduced from the data in Fig. 4,a are shown 
in Fig. 4,b. The manatee in which RTRM ceases to occur as the intensity of 
H., is increased will be observed in these figures. The disappearance of RTRM 
by the increase of H,, was interpreted to be caused by the fact that the effective 
magnetic field becomes positive when the intensity of H,, reaches some 200e. 
This value is remarkably smaller than that which was estimated in the previous 
paper (See Fig. 2 in [2-1], where the fact that the reversal of the magnetization of B was 
still occuring under H,, with an intensity as high as 40Oe. in partial TRM J; ,, was 
regarded as evidence indicating that H,>400¢e.). At that stage of study, the fixing 
of magnetization of B was supposed to complete during the time when H., is acting, 
namely at a temperature between 300°C and 250°C. But the data shown in Fig. 4,a 
suggest that such an assumption may not be correct: i.e. fixing of magnetization 
of B mainly takes place below 200°C. This relation is well illustrated in Fig. 5 
which shows the partial TRM charac- 
teristics of A’B-grains in strong fields. sail oe 
The partial TRM J +00, H,, (2 0» 
T>250°C is reversed even when H., 
amounts to 10000e. while the partial 
TRM J a %; uz, (Lo), T<200°C is nor- 
mal. This figure suggests that H, 
which magnetizes B in the reverse 
direction in Ji, y, » T>250°C, is due 
to the fixed part of magnetization of 
A’, and that the domain-fixing of B 
takes place after H,, is suppressed. 

In concluding this section, facts 
revealed at the present stage may be 
summarized as follows: (i) the 
magnetic interaction between A’ and 
B, responsible for RTRM, takes place -12 
within each grain of which the main Fig. 5. 
constituent is B: (ii) A, in the 


Partial TRM, J 7 +509 Hep (T,), of A’B-grains 
in strong magnetic fields. 
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separated A-grains and in the state of coarse attachment to B in AB-grains, has 
no essential meaning with respect to RTRM, and should be excluded from the 
study on its mechanism: (iii) A’ in A’B-grains is in a state of minute exsolution 
and its mean Curie point should be estimated to be about 320°C: (iv) the effective 
Hz working on B in A’B-grains at the time of the fixing of magnetization of B is 
of the order of 100e¢. 
§ 3. Magnetic Interaction in a System Composed of Two Ferromagnetic Materials 
with Different Curie Points 

It can be shown that the presence of two ferromagnetic constituents with 
different Curie points does not always produce a magnetically anti-parallel coupling 
between them. Therefore, it is presumed that the mutual relation of the two con- 
stituents in the case of RTRM must be of some particular configuration favourable 
for an anti-parallel coupling. 

In general, the total magnetization of a two-constituent system at a tempera- 
_ ture T, J(T), can be expressed as 


MT )=¢ea-JaT)+¢5-J AT), (1) 


where c is the volume concentration, and the subscripts A and B refer to the 
quantities of the two constituents A and B. Let us consider what will happen 
when this system undergoes a cooling from a temperature higher than the Curie 
point of A, 9.4, assuming 9.,>6.,, in a weak magnetic field kept constant 
throughout the whole process of cooling. Assuming that in a weak magnetic field 
the magnetization is proportional to H,, and J,, J(T) at 9,4>T>9,, is expressed 
as follows, 


MT) =Ca-JAT) = Ca AT) Jia T)* He (2) 


where 14,(7) denotes the constant of proportionality. For the sake of simplicity, 
JT) will be considered as fixed in the process of further cooling. Then, at T<9.,, 
if there is interaction between A and B, J(T) will be, 


HT) =carddT) S(T) Het enh T) Sil T (bez N-JAT)}, (3) 
since the intensity of the interaction can be presumed to be proportional to J,(T), 


where WN is the constant of proportionality determined by the nature of the spatial 


relation between A and B. 
When N>0, the interaction is negative and the magnetization of B will be 


reversed if the condition 


Hizx<N-JAT) . (4) 


is fulfilled; whereas when N<0, the interaction is positive, and therefore the 
effective magnetic field on B is more intense than H.;. Then the present problem 
of magnetic interaction reduces itself to the evaluation of N for a given configura- 


tion of the two constituents. 
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Let us consider some simple ideal 
cases. 

Case 1. Concentric spheres: when A is 
included in a concentric sphere of B, 
(Fig. 6) and J, is produced in the direction 
of H., that is oriented to the x-axis, the 
x-component of H, due to J, at (7, 9), 
provided that J, is uniform, can be 
expressed by 


20 —sin?0 
Hes =JarBar gs 9 5) 
where v, denotes the volume of A. The 
average value of H,;. over the volume of 


B, vz, is given by 


= ‘| Hy.dv. (6) 
Fig. 6. Model of co-existence of A- and B- Vp JrsB 
constituents in concentric spheres. 2 " 

It can be readily shown that the integral 


of the right-hand vanishes. There is, in such a case, no mutual magnetic interac- 
tion of significance for the present problem. 

Cases 2. Uniform mixing of A and B units: when a large number of minute 
units of A and B are mixed at random as in Fig. 7, the magnetic interaction may 
be approximately treated as being similar with the case of local electric field in 
dielectrics. In the course of field cooling, the effective magnetic field, Hz, which 
acts on B at the temperature where B is magnetized, can be divided into four 
components as 
Aug =H.2+Hi+Az+Hs, (7) ————-»> H ex 
where H, is the H, resulting 
from the magnetic poles pro- 
duced by J, on the outer 
surface of the specimen, H. is 
the field by the magnetic poles 
on the inner surface of the 
fictive spherical cavity and H; 
is the sum of the magnetic 
field by all the units of A 
contained in the cavity. (See 
for example, p. 92 [14}.) 
Though the units in the 
present case are utterly dif- 
ferent from those in dielect- 


oe ee Fig. 7. Model of co-existence of A and B-units i 
rics, it may be presumed that of random mixing. aT eaten 
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the calculation of H, may be treated in parallel manner with the case of the latter, 
provided that each A-unit can be approximated by a magnetic dipole oriented in 
the direction of H.. and the number of the units is sufficiently large. Then H. 
can be expressed by 


Hg =H.,—(D~32)cu-Jo (8) 


where D is the demagnetizing factor of the whole specimen in the direction of H.,. 
When D is small, as in the case of a long rod specimen, it is evident that H. 
becomes, 


4 
Hi.g=H.,+ gt Cala> } - Seat 


Therefore, in such a case, it may be expected that the interaction will be positive 

and that J» will be produced in the direction of H., even in partial TRM J7*x » 
where 7,>9_.,>T->9_,: namely even when H_,=0 at the temperature where the 

‘Magnetization of B appears. An experiment by the author showed that in a 

mixture of the grains of magnetite and nickel the latter grains are magnetically 

coupled by positive interaction with the former ones. 

In Néel’s model [4], the configuration is con- 

sidered to be composed of ellipsoidal B-grains a Mente B 

embedded in a spherical specimen containing the 
A-grains homegeneously. In his model, H.7 on B is 
expressed by, (Fig. 8), 


Hy =H..—( 37D) cu-Ja 


where D is the demagnetizing factor of the 
ellipsoid of B in the direction of H.,. When the 
specimen is not spherical but of a long-rod shape, wasn “anew anak ‘ 
as is in usual experimental conditions, the interac- of A- and B-constituents. 
tion will become positive in this case, too. 

Case 3. Parallel strips model: as seen in the above examples, it is by no 
means self-evident that the co-existence of two different magnetic constituents can 
results in the production of self-reversal of TRM. Although many other kinds of 
configurations of A and B can be considered, the evaluation of Hz in the case of 
parallel strips with alternate constiutents will be especially considered here, since 
the minute assemblages of minerals in rocks frequently take the form of inter- 
growth in parallel strips, like the Widmannstatten structure. In fact, the electron- 
microscopic ovservations described in the next section offered some evidence 
supporting the hypothesis that the present case is of a parallel strips structure. 

The model we consider is of parallel strips with alternate constituents as 
shown in Fig. 9: the plane of the strips is in the (yz}plane and the z-axis is in 
the upward direction vertically. The width of the strips and the height of the 
system are denoted by d and L. When such a system is cooled from a temperature 
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higher than 9.,4 in a magnetic ficld, H.,. 
applied in the z-direction, to a tempcra- 
ture T, where 9.,>T>9,,, the A-strips 
are magnetically fixed in the z-direction. 
_---Y For the present, let us assume for the 
- convenience of calculation that the 
Saeed oT Ree =e: magnetization of all the A-strips is fixed 
in the z-direction. The direction of the 


N 


magnetization of the B-strips, in the 
process of further cooling down to 
T<@.,n, is determined by the local Hig 
which is the sum of H,., and H, due to 
JT). Let us consider the cases, where 
Jex is oriented parallel (p-case) and 
where J, is oriented antiparallel (a-case) 
Fig. 9. nite etiam ins a. en mi to J,4, and compare the magnetic energy 
alternate constituents. in both cases. The necessary condition 

for the reversal of Jz is reduced to, 


E,< Ep, (9) 


where E stands for the total magnetic energy and the subscripts a and p designate 
the a-case and the p-case as defined above. Generally, the total magnetic energy 
of such a system can be represented to be 


E=Ewr+Enen, 
Evevw=EgwertEpveurtE ar, (10) 
Eren=E4,neu +Epnen , 


where Ever, Ena and E,, stand for the self energy, the interaction energy with 


H,, and the mutual interaction energy between J, and J», respectively. Enew and 
Exe are calculated by 


Erew = {TH dx dy dz, (11) 
1 
Exe=—} | P02 dF a (12) 


Now, since J, is assumed to be always fixed in the z-direction, Een can be dropped 


out of calculation. From Eq. (11), Ex.en for the unit area of the (xy)-plane is given 
by 


Exfew = —H..| ad pai et -L, (13) 


where the signs (—) and (+) refer to the a@-case and the p-case. The magnetic 
self energy of the system Eu... is calculated approximately by a similar method as 
shown by Kittel [15]. In the a-case, the 2-component of the magnetic field just 
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bellow the surface due to the free poles 
on the surface is illustrated by Fig. 10. 
Expressing this distribution by the 
Fourier expansion, the approximate 
solution for the Laplace equation is given 
by, provided L>d,+-d,, 


H.=3i\a,+cos nkx-e"*, (14) 
0 
where k= 2x and a,, are the Fourier 
datd, a‘ 
coefficients namely, 
2n 
a= J 4s LJn-dy—Jardal, 
¥ in an 
Fig. 10. Distribution of magnetic field just mre (Jat+Jn)sin dyt+d,~ (15) 
bellow the surface in parallel strip d . 
model. (a-case.) Putting (14) and (15) into (12), we can 
obtain Eur as 
rl. .* . *s " 1/ .. nad Pr iag 
Ever = (datds)? {JardacZoJnedy)* + A(Jac= Ju) n (sin 704) [ie hd 
dp/2 d4/2+dR/2 
x aay & > : +<cosnkx > ae ye v 
| 0 dB/2 J 
; ~P_ 1 datdy ... nndx 
<cosnkx > Rage Path ad HP ae 
7 tareanls 1 datdn i. mnds 
<cos nkx Be os alet aiie POo (16) 


where, since L>d,+dp,, the integration with regard to z is carried out from —« 
to 0 for all the terms except the first, and doubled to account for the effect from 
both of the top and bottom surfaces. That the signs in parenthesis in the above 
formula refer to the p-case can be shown in a similar manner. The final result 
for the condition of the reversal of J, from Eqs. (9), (10), (13) and (16), becomes 


ae: Cs a ee ony 
Hee<2Jay 46 Ant 1 +e |’ 
iL _ dz 17 
where A= ads” and c= i (17) 


This relationship is illustrated in Fig. 11. It is assumed in the present model, as 
has already been remarked, that the magnetization of all the A-strips is fixed in 
the direction of H,,. However, unless the A-strips are magnetized to saturation, 
such an assumption is not warranted. Therefore, in order to utilize the condition 
(17) for the actual case, some statistical consideration on the arrangement of the 
magnetic moment of the A-strips will be necessary. It will be unnecessary to 
account for the case when H., is applied in the direction perpendicular me oe 
(yz)-plane, because it can be expected that the thin strips will hardly be magnetized 


in the x-direction (Fig. 9). 
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° 
100 C= ery, 10 i tho Vico 


Fig. 11. Plot of Eq. (17). Bellow the curves the anti-parallel 
arrangement of J, and J, is energetically more favourable; 
above the curves the parallel arrangement is more 
favourable. 


An experimental study on the model somewhat similar to the present case 
has recently been reported by Grabovsky and Pushkov [16]. They showed, by 
direct measurement of the local magnetic field, that the relation H,~-H,, can be 
fulfilled in the limited region near the plate-like specimens of magnetite when 
they acquired TRM in the earth’s magnetic field. The dimension of their specimens 
were of the order of 10cm. 

§ 4. Electron-Microscopic Observations on the Intergrowth of the Two 

Constituents, and Estimation of the Phase of the A’-constituent. 

According to the calculations described in the preceding section, it may be 
expected that the parallel strips model can be responsible for the self-reversal of 
TRM under appropriate conditions. It will be reported in this section that such 
a structure was actually observed in the RTRM-producing grains (A’B-grains) of 
the Haruna rock with the aid of the electron-microscope. Figs. 12, a. and b. are 
the electron-micrographs of etched surfaces of A’B-grains in the Haruna rock. 
The observations were carried out by the methyl-metacryl alminium replica 
method. The strips seen in the figures were produced on the polished surfaces 
by etching with HCl for about 15 minutes or with HF for 1~2 minutes. The 
width of the strips is approximately 0.2, while the length ranges apparently 
between ly and 5,. This scale may be fine enough to be comparable with that 
of one or several elementary magnetic domains. It seems likely that these strips 
are more or less chained to one another along specified directions to form fairly 
long and narrow strips. We may, therefore, presume that this structure consists 
of the intergrowth of the two constituents which makes these grains acquire the 
RTRM characteristics. 

As for the phase of these strips, so long as we assume them to be A’, there 
are, as has been stated in § 2., two possibilities of identification : namely the Ti-rich 
Ti-Mt and the II-Ht solid solution. Judging from the known fact that reagents like 
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Fig. 12. Electron micrographs of polished surfaces of A’B-grain of Haruna ferromagnetic 
mineral, after etching with HCl for 15 min., showing intergrowth of A’ in echelon 
into B. 


HCl and HF etch the Ti-Mt more effectively than the II-Ht phase, it was at first 
inferred that the phase of the strips should be the Ti-rich Ti-Mt [17]. But such 
an approach may be not so reliable as other methods like the crystallographic 
one. A powder method X-ray analysis, with the aid of NORELCO spectrometer, of 
A’B-grains showed no trace of the Ti-Mt phase. However, an X-ray analysis on 
a single crystal of A’B-grain showed that, in the oscillation photograph, the Debye- 
Scherrer rings of the II-Ht phase could also be observed, together with the intense 
spots by the single crystal. Fig. 13 is an example of the oscillation photographs 
taken by using Mo K-radiation. The Debye-Scherrer rings, though faint, have 
non-uniform intensity. This fact may 
indicate that, in a single crystal of A’B- 
grain, some parts of the II-Ht phase 
exist in a state that can cause the 
powder-like diffraction with a certain 
degree of orientation. Although the A’- 
constituent originating the RTRM, the 
strips observed by the electron-micro- 
scope, and the constituent originating the 


Fig. 13. Cscillation photograph of A’B-grain Debye-Scherrer rings in the single crystal 
of Haruna ferromagnetic mineral, - 3 
using unfiltered Mo K-radiation, have not been completely masa to be 
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the various data obtained for the A’B-grains, to consider these to be the same. 
Therefore it will be concluded that the phase of the A’-constituent is Il-Ht, with 
the Curie point higher than that of B. 
§ 5. Application of the Parallel Strip Model to the Self-Reversal of TRM of the 
Haruna-Type 

A theoretical consideration of the mechanism of RTRM of the Haruna-type 
will be discussed, in this section, in terms of the parallel strip model treated in § 3. 
As it was verified that the ferromagnetic grains responsible for the production of 
RTRM of the present case are the A’B-grains, it would be sufficient to explain the 
RTRM characteristics which are summarized in Figs. 4,a. and 4,b. The rest of 
the ferromagnetic minerals present in the original Haruna rock, namely the A- 
constituent, may be considered as being non-essential with respect to the RTRM 
phenomenon. The principal features of the hereby proposed model are as follows. 

The host mineral in the A’B-grains is B, a ferromagnetic II-Ht solid solution 
whose Curie point 9., is about 230°C, and there are fine strips of A’ whose Curie 
point 9.,’ ranging around 320°C. In the process of cooling in a weak magnetic 
field, H,,, the TRM of the A’ strips is produced at a certain temperature just 
bellow 9,4’, so that there would appear the parallel-strips effect which was examined 
in § 3. In other words, the magnetization of the A’-strips would exert a reversely 
oriented magnetic field H,; upon the part of B which is between the A’-strips. 
The parallel-strips effect on the part of B 
which is not between the A’-strips may be, 
in the first approximation, considered to be 
weak enough to be neglected. In order to 
distinguish these two parts of B, the latter 
will be called B’ (Fig. 14). If H, is stronger 
than H,,, the TRM of B will be produced 
in the reverse direction when the A’B-grains 
are cooled through 9,,, and if this reverse 
TRM of B is stronger than the normal 
TRM of A’ and B’ at the room tempe- 
rature, the TRM of the whole A’B-grains is 
reversed. 


Fig. 14. Structure of A’B-grain. 


It may be possible that the B-constituent which is between the A’-strips has 
the properties of fine grains as the A’-strips, since it is finely divided by the A’- 
strips. As for the nature of these two parts of the B-constituent, therefore, we 
may presume the following characteristics: their structure insensitive properties 
such as the Curie point and the spontaneous magnetization are identical but their 


structure sensitive properties such as the coercive force, the magnetic susceptibility 


and the TRM characteristics are different. Then basic assumptions on the nature 


of each constituent in the present discussion will be: 


i) A’ and B have the typical TRM characteristics, such as the large Q-value 
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and the saturation phenomenon of TRM in the magnetic field as weak as the order 
of 100O¢. [18], [19]. Because the size of these constituents may be as sufficiently 
small as to become the case of the single-domain model of TRM by Néel [20], which 
could very well explain the typical TRM characteristics such as mentioned above. 

ii) B’ has not the typical TRM characteristics, because the size of this 
constituent is fairly large. Its Q-value is small, say less than unity, and there is 
no saturation phenomenon in TRM within the magnetic field of the order of 100e. 

At a temperature T, the fixed part of the magnetization of the A’B-grains 
can be written as 


Fixe T)=0.4:+J4',x0a(T )+¢2*Jp, axea(T) +n] x’ 2x0a(T), (18) 
Javiixe T)= 8, 4'*-Jea(T)* Hi, 4 

Ix, xe T)= 8, 2-JsT)-{H..—Ha}, 

Tr’. xe T)= 8, n'*Jen(T)* Hers 


where H; is the mean value of H, due to Jy and s, denotes the domain fixing 
ratio defined generally by, 
Taxed T)= 8,°J, T)-H. (19) 


Now, the direction of the plane of the strips are supposed to be equally distributed 
in the (xy), (yz) and (zx)}planes (Fig. 15). z 

As stated in § 3, when H_, is applied in the 
z-direction, the A’-strips and B lying in the 
(xy)-plane do not contribute to the magnet- 
ization. Therefore, the amounts of A’ and 
B relevant for the magnetization are 
reduced to 2/3-c, and 2/3-c, respectively. 
As also mentioned at the end of § 3, in order 
to estimate the actual value of H, by the 
parallel strip medel, the effects due to the fact 

that not all of the A’-strips are magnetically Fig. 15. Equal distribution of parallel strips 

in (xy), (yz) and (zx)-planes. 


fixed in the direction of H,, must be taken 
into account. Since the concentration of the A’-strips lying in the (yz) and (z2x)- 


planes is presumed as being 2/3-c,’, the fraction of the A’-strips of eg magnetic 
moment is fixed in the direction of H,, may be estimated to be 2/3. os —, where 
Jw is the intensity of magnetization of A’ as a whole fixed in the direction ar 77... 
The magnetic moment of the remaining fraction ( 1-4 /3- a =) of the A’-strips 
should be considered to be equally distributed in z2- and —z-directions. Then, in 
the first approximation, the amount of B which is under the influence of the 
parallel-strips effet may be regarded as reduced to ae xX 9 at = , and the rest of 


B, namely 3ea( 1- z whip 1 -) may be regarded to be under the effect of H,, only. 


Since the average value of the ratio d;/d,’ is approximately estimated to be five 
from Fig. 12, the intensity of H, in the present case may be estimated as being 
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nearly equal to J,4’, according to the curve in Fig. 11, provided that the magnetic 
moment of the A’ is of the order of spontaneous magnetization. 
Then, the average value of the effective field, H.g, over all the B-constituent 
which can contribute to the magnetization of the A’B-grain may be obtained as 
Hg (Hus—Ioa) 9 4 + Hee'( 1-978), )= Hee 9 Ie (20) 
The reversible part of the magnetization of the A’B-grains at the room 
temperature 7, is, on the other hand, expressed by 
Jreversnie( To) = An T.)*€a'*Jeh T0)* Hees (21) 
provided that the contribution from A’ and B can be neglected by the assumption 
i). 
From Eqs. (18), (20) and (21), the Q-value of the A’B-grains at the room 
temperature 7, will be written as 
Ce Sun Isl T.)) Hea 3 Ju(To)| +¢4'*H, Sy a'*Jsa Ty) +€n'* Hee Scn'* Jen To) 
Q(HA,.)= ; Crh’ To) Isa To): Hee 
where J4(T)2)=Sp4':J:4(Ti2):H.. and T,, denotes the domain fixing temperature 
of B. Although it is expected that J,,4(To) is larger than J,47T >), since 9.4°>90.2 
[21], the difference in J, expected from the difference in the Curie points may be 
small. Therefore, we may put, in the present approximation, J,4(T)=J.«T.) in 
Eq. (22). We may also put s,”’~A,(T))/2 by the assumption ii). The value of A4;(T) 
can be obtained as wax 800 by the definition of 4 (Eq. (2) rewritten for B’.) and 
the data in Fig. 4,b. The value of J,4(7,s) will be about 40 emu/cc, provided 
that A’ is a ferromagnetic I]-Ht of which J,/T) curve is known to be linear, and 
J,a(T.) is nearly equal to J,,(T)) which has been determined to be about 120 
emu/cc previously [3]. 
Putting the above relations and values, together with another relation 
5¢4'=Cz (Fig. 12), into Eq. (22), we can obtain 


(22) 


Q(H,,)=4 x 10%+¢4/*S, x(1—60- 8, 4’) +8x 10? +€4/*S,4'+1/2. (23) 


In the next place, the saturation effect of the TRM of A’ and B (assump- 
tion i)) will be expressed approximately by 


Insel #1) = M(1—e-#/*), ; (24) 


where M and hare the saturation value of TRM divided by J, and the intensity of 
the magnetic field in which the TRM attains the (1—1/e) of its saturation value. 
From Eqs. (19) and (24), s,4’ and s, will be expressed by 


My: 
8, ¢= 1—en~ex/har 
H,, § ) 


M 3 
Spz= | Hag | (l—e- !4egV/“8), (25) 


where the absolute value must be used for H.g Since it can be negative. Then, 
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the final expression for the Q-value becomes 


LA 1—e-| Meg fir i er : 
QWHe.)=AX 10% C47-My "9 F He 8X10 -c4r Myr! A “chy 9, 
where Hep =H.—60>Mg(1—e-"er/ ha’), (26) 


In Eg. (26), the value of ce,’ may be put as 0.1 from Fig. 3,a. The values of the 
parameters My’, Mx, hy and h,; are inaccessible at present. The saturation value 
of TRM was experimentally studied by Melle. Roquet [22]. According to the result 
of her study, M value may be estimated to be about 1/8 for magnetite. As for 
the value of h, Nagata [19] and Thellier 
[18] have reported that some natural 


5 Q-value 


magnetite shows the non-linearity 9 
between the amount of TRM and the 

intensity of H., when H,, is of the order 

* Of ou M00. Since the ferromagnetic oy 
materials in the present case are the 
ferromagnetic II-Ht solid solutions, there 
seems to be no experimental data 
available for the estimation of these 
parameters. The Q(H.,) curve in Fig. 16 
has been calculated from Eq. (26) by 
putting, for instance, -30 


My=1/3, hy =2, 
M;=1/5, hz=10. 


In Fig. 16, the experimental Q(H.,) values 
for the A’B-grains of Haruna sample, Fig. 16. Field dependence of Q-value, calculated 
obtained from the data in Fig. 4,b are al bey Pee oe (thea 
also plotted for comparison. experimental Q-values. A’B-grains. 
§ 6. Other Rocks Having the Similar Composition of Ferromagnetic Material as the 
Haruna Pumice 

There have been found, so far, 10 other rocks which contain the ferromagne- 
tic Il-Ht solid solution together with Ti-Mt as in the Haruna pumice. The 
summary of the general characteristics of some of these examples reported in the 
recent paper by Akimoto [21] is reproduced in Table 1 with the permission of 
the author. Some tentative results of the examinations conducted on these samples 
with regard to the magnetic interaction in the process of the development of TRM 
will be described in this section. The examinations consisted of the following 
procedures, which had already been applied to the Haruna sample. 

1) The thermal change of susceptibility «(T) and the modes of development 
of the total TRM in weak magnetic field and its decay by heating in a non-magnetic 
space. By a study of these processes, the fundamental characteristics of the sample 


with regard to TRM may by clarified to some extent. 
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2) The mode of the development of the partial TRM To, 4m, & T),. where 
T,>9,,. This process was conducted in order to determine the sense of the 
magnetic interaction between A and B, if it exists. In this process, no magnetic 
field but H, due to J, works on B at the domain-fixing temperature of B. 

3) Thermo-magnetic separatior : by this operation, it was tested whether or 
not grains mainly consisting of the B-phase as in the Haruna rock could be found. 

4) The «(T) and the modes of the development and the thermal decay of 
total TRM of the thermo-magnetically separated B-grains, which correspond to the 
A’B-grains in the Haruna sample. The magnetic interaction within the B-grains, 
if such exists, should be observed, in this process more effectively than in 1), because 
the effect of the non-essential A-constituent is eliminated. (See, for example, 
Fig. 2.) The TRM characteristics, such as the Q-value, of the ferromagnetic Il-Ht 
phase are also obtained by this process. 

5) Electron-microscopic observations on the polished surfaces of the thermo- 
magnetically separated B-grains. In these observations the Triahol-alminium-replica 
method was applied on the surfaces of the samples etched with HF for 1 minute. 

The results of these examinations obtained on each sample is summarized 
bellow. For reference, the case of the Haruna rock, as has already been stated, 
will be summarized first as follows. 


Table I. Properties of Ferromagnetic Ilmenite-Hematite Solid Solutions. [21] 


|—— = a aft 


: . | LE 9 = ay * ef ee 4 
|_ Locality , Rock vars Feo Tio, ie Lf sigromenas 6. . I(T) 
) dacite | . oe . .| 
Haruna | Sumice 26.9 33.2 39.9 5.480A 55°08’ | 250° il 
er |” dacite ae = ea ey verte ay 
BHO | Sitetetone 1 | 5.483A 55°02’ | 230° |~21 ,, 
Sokota iron sand | 5.483A 55°01’ 190° ; 
Dean tia ~ hyp. horn. , +r = 
Towada | Gyr Ornnice 19.63 38.37 42.00 | 5.491A 54°59’ | 100° | 24 ., 
Okub: ‘ HS et eee F | 
eee a mic dacite obsidian 20.34 40.88 38.78 5.495A 54°58’ | 10°) 21, | 
Asama dacite pumice | 16.97 40.66 42.37 5.498A 54°56’ | 100° | 18 ., | 
j i horn. “mica : “~ > ¢ ~ i Pearse) ; ig a 
Himesima| ye nesite | 20.70 35.53 43.77 5.489A 54°56’ | 190° | 29.» | 


* in mol. % 
Haruna Rock : 

1) The magnetization of B is produced in a reverse direction, during the 
cooling process in a weak magnetic field, and the total TRM is also reversed at the 
room temperature. (See Fig. 20, [1-4].) In the thermal decay curve of the TRM, the 
anomalous increase of TRM was observed at about 450°C. (p. 161, [Rock-Magnetism]}) 

2) The modes of the various partial TRM are shown in Figs. (11), (12), (13), 
(14), (15), (16), (17) and (18) in [1-4}. It is quite evident that the magnetization of 
B is reversed in some of the partial TRM. 

3) The A’B- 


grains were separated out by the thermo-magnetic separation 
[2], [3]. 
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4) The mode of the development of the total TRM for the A’B-grains is 
shown in Fig. 2. The reverse TRM of the B-constituent was about 30 times more 
intense than that of the original sample. 

5) The well developed A’-strips were observed. (Fig. 12) 

Asio Rock : 

No examination has been newly made, since the sample was used up in 
the experiments reported in the previous paper [3]. The locality where this 
rock is exposed is said to have sunk under an aartificial lake constructed for 
power development. The results reported in the previous paper are summarized 
as: 

1) The total TRM of the original sample was normal, but it was verified 
through a study of TRM of the single grains that the magnetization of B is reversed 
within the total TRM. 

2) The partiai TRM /;"" Mn, ( T,), H..<3.5 Oe. was proved to be reversed. 


3) The results of thermo-magnetic separation is quite similar to that of the 


Ix 10° emu Z, 


Haruna case [3-2]. Most of the grains with 
the thermo-magnetic separation temper- 
ature below 250°C showed perfect RTRM 
characteristics. 

Sokota Iron Sand: 

The self-reversal of TRM of the 
Asio-type of the iron sand at the Sokota- 
mines was, as mentioned already, dis- 
covered by Saito [13]. 

1) The experiments for the total 
TRM were carried out by the above 
author, according to whom the reversal 
of TRM of B, which is not large enough 
to overcome the normal TRM of A was 
verified. The anomalous increase at 
about 500°C in the heating curve was also 
observed. 

2) The modes of the various parti- 
al TRM were also studied by the above 
author and it was proved that the 
characteristics are quite same as the 
Haruna and Asio cases. 

3) The thermo-magnetic separa- 


Fig. 17. Characteristics of thermo-magnetically 


tion conducted by the present author at ee iad Ss eit OF Sokdta trod badd, 
250°C, on the sample obtained by in H,,=2.0 Oe. 

. be Full circles: «(T)xX He, curve, 
courtesy of Mr. Saito, showed that Hollow circles: development of total TRM. 


B-grains, or A’B-grains are present. Broken line: thermal decay of total TRM. 
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4) The mode of the development of the total TRM of the thermo-magnetically 
separated B-grains shows that they are mainly responsible for the RTRM character- 


istics of the Sokota iron sand. (Fig. 17) 
5) In Fig. 22,a the well developed strips can be observed. 


Towada Rock: 

1) As seen in Fig. 18,a, the TRM of B is produced in the normal direction 
by the process of total field cooling. This type of mode of development of TRM 
will be called the Towada-type [10]. The anomalous increase in the heating curve 


was observed at about 430°C. 

2) In the process of the development of the partial TRM 1 a. OA be shown 
in Fig. 18,b, the intensity of the remanent magnetization was observed to decrease 
as the specimen was cooled through 9,,. The results of 1) and 2) may be interpreted 
as showing that the magnetic coupling in the Towada sample works anti-parallelly, 
but that it is so weak that its effect does not appear in the process of totai field-cooling. 

3) The presence of B-grains was affirmed by the thermo-magnetic separation 
at about 150°C. 

4) The experiments on the B-grains are shown in Fig. 18,c, in which the 
normally oriented magnetization of B will be seen. 

5) Very scarce traces of the strips were observed. (Fig. 22,b) Therefore, 
the weak interaction assumed above might be explained by the effect of the rarely 
found strips. 
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Fig. 18, a. Fig. 18, b. Fig. 18, c. 


Steinar he one Development of partial TRM Characteristics of thermo- 
inera ; i 
s of Towada rock in Jii0.8.00e(T) of Towada ferro- ™agnetically separated B- 


Fi oi= 2.006, : 
Full circles: «(T)xH,, curve, | ™agnetic minerals. grains of Towada _ ferro- 
Hollow circles: development of magnetic minerals, in H,, 
total TRM. iret 
Broken line: thermal decay of Full circles: «(T)x H,, curve. 
total TRM. Hollow circles: development 
of total TRM. 
Broken line: thermal decay 
of total TRM. 


Okubosawa Rock : 
1) The mode is the Towada-type. (Fig. 19,a) 
2) The mode for the partial ‘CRM J®™  (T) is shown in Fig. 19,b. The 


500, 8.006, 
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remanent magnetization increases as the temperature decreases. It was not possible 
to determine whether this increase in TRM is due to the production of the 
normally oriented J; by the positive interaction with J,, or due only to the increase 
of the spontaneous magetization of A. 
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Fig. 19, 


a. Characteristics of ferromagnetic 
minerals of Okubosawa rock in H,, 


= 4.008. 
Full circles: «*(T)xH,, curve. 


Hollow circles: development of total 


Broken line: thermal decay of total 


TRM. 


b. Development of partial TRM, 
Ji) s.00e(T), of Okubosawa ferromagne- 


tic minerals. 


5) The comparatively rare occurence of 
obscure strips was observed by electron- 
microscopic observations. (Fig. 22, c) 
Asama Rock : 

1) Although the intensity of magneti- 
zation is weak, the comparatively intense 
normal TRM of B will be seen in Fig. 20,a. 
This mode may be classified with the 
Towada-type. 


Ix10°” emus, 


Fig. 20. 
a. Characteristics of ferromagnetic minerals 
of Asama rock in H,,=4.00e. 
Full circles: “(T)xH,, curve. 
Hollow circles: developmedt of total TRM. 
Broken line: thermal decay of total TRM. 
b. Development of partial TRM J}i0, s.o0e(T ) 


of Asama ferromagnetic minerals. 


2) No effect of the magnetic interaction was observed (Fig. ,20,a) in this 


case, since the partial TRM /’” 


(T) showed only a very slight increase below 


500, 8.00¢. 


150°C. 


3) B-grains were separated out by the thermo-magnetic separation at about 


150°C: more than 90% consisted of the B-grains. It will be of a petrological 
interest that the major constituent of this sample is [l-Ht, instead of Ti-Ml. 
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4) No distinct strip-shaped intergrowth was observed. In Fig. 22,d, some 


etched spots will be observed. 


Himesima Rock : 
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a b 
Fig. 21. 
a. Characteristics cf ferromagnetic minerals of Himesima rock in H,,=4.0 Oe. 
Full circles: «(T)x#H,, curve. Hollow circles: development of total TRM. 


Broken line: thermal decay cf total TRM. 
b. Development of partial TRM J?” ra (T) cf Himesima ferromagnetic minerals. 
OO, 4. et. 


Characteristics of thermo-magnetically separated B-grains of Himesima ferromagnetic 
minerals in H,,=4.0 Oe. 

Full circles: «(T)xH,, curve. 

Hollow circles: development of tctal TRM. 

Broken line: thermal decay of total TRM. 


d. Development and thermal decay of partial TRM, /*!” 


21 
0 ; 350, 8.00¢. 
separated B-grairs of Himesima ferromagnetic minerals. 


(1) of thermo-magretically 


Fig. 22. a. Fig. 22,b. 

Electron micrograph cf polished surface Electron micrograph of polished surface 
of thermo-magnetically separated B-grain of thermo-magnetically separated B-grain 
of Sokcta iron sand, after etching with HF cf Towada ferromagnetic minerals, after 
for 1 minute, showing well developed etching with HF for 1 minute showin 
exsolution lamellae, scarce exsolution lamellae, ; 
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Fig. 22,..¢. Fig. 22, d. 

Electron micrograph of polished surface Electron micrograph of polished surface 
of thermo-magnetically separated B-grain of thermo-magnetizally separated B-grain 
of Okubosawa ferromagnetic minerals, of Asama ferromagnetic minerals, after 
after etching with HF for 1 minute, etching with HF for 1 minute, showing 
showing obscure exsolution lamellae. undeveloped exsolutions. 


Fig. 22, e. Fig. 22, 1. 

Electron micrograph cf pclished surface Electron micrograph of polished surface 
of thermo-magnetically separated B-grain of ferromagnetic grain of Adirondack rock 
of Himesima ferromagnetic minerals, after No. 21, after etching with HF for 1 
etching with HF for 1 minute, showing minute, showing well developed exsolution 
well developed exsolution lamellae. lamellae. 


1) The mode obviously belongs to the Towada-type. (Fig. 21,a) 
2) The partial TRM, J*” — (T), increases as the sample is cooled after 


vo 4.0 Oc. 
the suppression of H.,. (Fig. 21,b) 
3) B-grains were separated out by thermo-magnetic separation at about 
250°C. 
4) The TRM of B was in the direction of H,,. (Fig. 21,c). For the 
B-grains of this sample, the mode of development of partial TRM Rey T), was 
also measured. (Fig. 21,d). In this process, an anomalous increase in the remanent 


magnetism after the suppression of H,, was observed bellow 100°C. Whether or not 
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this increase is caused by the positive magnetic interaction between the constitucnt 
with lower Curie point and the constituent which acquired TRM before the sup- 
pression of H,,, is not clear at the present stage. However, the possibility of the 
occurence of such an effect seems to well deserve further examination. 

5) It may be of interest that the patterns of the strips observed (Fig. 22,e) 
are similar with those in the cases of RTRM. In order to clarify the mechanism 
of the magnetic interaction in this case. it would be desirable to study the nature 
of the strips; their phase, their Curie point, etc. 

Finally, the ferromagnetic mineral-grains extracted from the Adirondack 
rock, which was offered by Dr. Balsley to Prof. Nagata, were examined by means 
of the electron-microscope. (Fig. 22,f) The rock-sample was No. 21 which was 
reported to have reverse natural remanent magnetism [11]. According to Balsley 
and Buddington, this reverse natural remanent magnetism is caused by the self- 
reversal of TRM, though the chemical composition of the participating ferromagne- 
tic constituents differs from that of the Haruna case. Judging from the chemical 
composition reported by the above authors and the results of the X-ray analysis 
conducted in our laboratory, the host in Fig. 22,f may be presumed to be a 
titanhematite and the strips to be another titanhematite containnig less Ti than 
the former. 

§ 7. Conclusions and the Further Problems 

When more than two ferromagnetic materials are present in rocks, it is 
generally expected that there exists some magnetic interaction between them, 
especially during the process of development of thermo-remanent magnetism. The 
characteristics, such as the sense and the intensity, of the interaction should 
naturally depend on their state of co-existence. It seems, therefore, quite probable 
that such presence generates various phenomena in thermo-remanent magnetism, 
which would otherwise hardly occur. 

In the Haruna rock, which is known to have the self-reversal characteristics 
of thermo-remanent magnetism, the titanomagnetites ard the ferromagnetic solid 
solutions between ilmenite and hematite are contained as co-existing ferro- 
magnetic ingredients. In the present paper, a detailed study was made, from the 
above mentioned viewpoint, on the actual mechanism of the magnetic interaction 
responsible for the particular character of the reverse thermo-remanent magnetism 
in this rock. The various data obtained by magnetic, crystallographic and 
electron-microscopic examinations on the A’B-grains, which were proved to be 
responsible for the phenomenon, demonstrated that the magnetic interaction should 
be assumed as being between two different ferromagnetic ilmenite-hematite solid 
solutions and that the titanomagnetites should be assumed as having little effect. 
A theoretical consideration also showed that the reverse thermo-remanent magnet- 
ism of the Haruna rock may be explained by the magnetic interaction between 
the two constituents intergrown in such a state as observed through the electron- 
microscope, provided that the fine strips of the ilmenite-hematite phases can acquire 
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a very strong thermo-remanent magnetism in weak fields. 

Examinations were conducted on several other rocks having similar assembl- 
ages of ferromagnetic minerals as in the Haruna rock. It was shown, through 
these examinations, that there actually exist various degrees of magnetic interac- 
tion in natural ferromagnetic substances. 

Though the studies in the present paper were confined to the rocks 
containing the ferromagnetic ilmenite-hematite phases, it would be_ highly 
desirable to study the nature of the magnetic interactions which take place 
between titanomagnetite phases, since it is known that the latter phases 
are of much greater importance than the former from the viewpoint of rock- 
magnetism as a whole [21]. In fact, several cases where magnetic interactions 
between titanomagnetite phases are expected to be playing an essential role 
have been reported. For example, Kawai et al. [23] have reported the occurence of 
an interesting phenomenon where the natural remanent-magnetism of some rocks 
changes its sign when heated in a non-magnetic space. The anomalous increase of 
thermo-remanent magnetism during its thermal decay process which was found in 
the ferromagnetic minerals from Niisima (p. 161. [Rock-Magnetism]) and in the A- 
constituent of the Haruna [1-4], Sokota and Towada rocks has also been inter- 
preted as caused by an anti-parallel coupling between titanomagnetites. In 
order to utilize the natural remanent magnetism of rocks for palaeomagnetic 
purposes, it was proposed [10] that the natural remanent magnetism should be 
proved as having the very normal characteristics of thermo-remanent magnetism 
and as being free from any anomalous phenomena like those mentioned above. 
Therefore a more thorough clarification of the nature of magnetic interactions 
in natural ferromagnetic substances will be of great significance from the viewpoint 
of palaeomagnetism. 
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The Thermo-Magnetic Properties and History of Some Plutonic 
Rocks from the Leinster Granite, Ireland 


By Horace MANLEY* and David J. BURDON** 
(Read Oct. 26, 1949) 


Summary 


The susceptibility and remanent magnetisation have been determined for 
seven representative samples of the granodiorite, granite and basic enclaves which 
form the bulk of the Arrigle Complex, lying at the south-west termination of the 
Leinster Granite. The determinations were made on the rocks in their natural 
state and again after they had been heated to 600°C and cooled in the earth’s 
magnetic field in London. From these have been calculated X, the ratio of the 
susceptibility after and before heating; S, the ratio of the thermo-remanent 
magnetisation to the remanent magnetisation ; and Z, the ratio of Q to T. 

For the granodiorites, X averages about 30 and S about 50; this indicates 
that new magnetic minerals were developed by heating and suggests that the 
granodiorites, as now constituted, have not cooled through their Curie point. For 
the granites, X averages 1.3 and S about 32; this indicates that no new magnetic 
minerals have been developed by heating, but that the remanent magnetisation of 
existing magnetic minerals has been greatly increased. This proves that the magne- 
tic minerals of the granite have formed below their Curie point; the petrography 
indicates that magnetite has formed or reformed along with other deuteritic changes, 
the most important of which is the change of hornblende to a biotite-quartz symplec- 
tite. Of the basic enclaves, the gabbroes show X equal to 1.6 and S equal to 2.7; 
this proves that the gabbroes have crystallised and cooled from above their Curie 
point in the normal manner of igneous rocks. Another enclave, supposedly a 
porphyretic dolerite, shows X equal to 3.9 and S equal to 21; either the dolerite 
has suffered much low-temperature metamorphism of its magnetic, or potentially 
magnetic, minerals, or it is a metamorphic rock closely resembling a dolerite. 

Thermo-magnetic investigations, c i i i 
outlined in this paper ol eae Pgh intr Pi a oe penvera th x: a 

n of the 
thermal history of many metamorphic-plutonic-igneous rocks. 
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1. Geological Setting 
The Caledonian granites of the Leinster chain form a continuous outcrop 
extending south-south-west for some 70 miles from Dublin to the village of Inistioge 
in Co. Kilkenny, Fig. 1, A. There is then a gap of some six miles between the end 
of the main granite outcrop and the present cover of unconformable Upper Old 
Red Sandstone on the west. Within this gap there are five outcrops of granitic 
rock. Four of these resemble closely the granite or granodiorite of the main 
Leinster chain; but the most westernly outcrop, in part covered by the Old Red 
Sandstone, consists of several different rock-types Fig. 1, B. It has been named 
‘The Arrigle Complex’ from the small river which flows over it; the rocks 
described in this paper all come from the Arrigle Complex. 
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Fig. 1. Generalised geological map of eastern Ireland, showing 
the relationship of the Arrigle Complex to the Leinster 
Granite. Details of the geology are shown on the second 


map. 


— 


ee ee ee ew ee nadie Ue ae Ey amy s alee} 
. : ’ 


voawol .| .d bas varwaM .H 


goitreqoig oiisagxM SO . 

; 2oitT9qGO1"d gees - ey 
he Thermo-Magnetic Properties and Higkeneseliome Plutonic 
: Rewks f) cay . & d@evadboa thro SrAand we . 
patio oitengsm-orrtedT &- 4 : oa 
veoretA2e?, oat, to. yrojeid Bre Boch a 
edqargoildie | bes ainomsgbe w +B ‘a 
nnd tan a . : 
21 ods lo esting tg nsinobsieD odT ‘= 


gorse evouniiaes ¢ mriet aisdo rotent iol 9 yy 
sgoiieinl to oystliv ony of ciidud mort eolien oy ‘Stno8 rot teow-ivoe-diuoe gaibastrs = 


bao ors néswted. elim! xie arabe Ie) G&R) B mae gadt fhe! ei, “GRETA? fhe 
DIO HeqUesidaniratao cur td. 2902 IOIRW, wat de AR, GOTAIUS. STB 
aiffitittg 20 Sqoroitfo: vit exe !owWdh gan. ait? AMIE. qtRAwW. Al 92. 29QRPGRE, B 


nist SY Focetitibéasrs!-t0r etiestanedd, loamy. eldmoene oad, te aud a 
hae HO Sitroed obeievcotinggh ad .qoun!9e:; sacar a 2 | 
ROMME Aeod ett 1 | ahU>ng!® eaqyiielees saat oMib,, leyeap, AM Sie? RAO UES 


| SSCS at advo awel hdéitw hewia.disme ad? mop} ox? ~ 


canpmetimtin to tie a snork ore@ Hs ae gi podixoneb | 
cha ane ~ cote coin BRSMOD ,. AOAAAK ae shout bet indicates 


le Reoagan SR oo VES ic coke potne 
ee gre <3" : 


einer sis heal he 


etofiucy *tagteetet 
wh mantras @ 
Cn See 


hw O4her Sw 


sTy 
Wie wae 
ge? . 
pease 
has 8 
eae , ee : =a Hic > s q 
XS i r 7 uy SNe. shits 
~~ _ ~ aoe 
Ti x oy Se he: ; ws: —_ pre ae 


> sce So 


arf agai Meer ote 
EASA << 


on . of 
2 2&8 SENG 


thermal bis 


AM 
- 


rie 


ee ~| e 


Te. LA gs. Gyr ra, 


“> eee 220A wrodwlt seed. HS eetiwle ba ccartsaow) atloxge thawed T ad T 


little fine-grained -erapaetisiod: reise ovilo. etzienos xviqmoD a 
magnetite oveurs ir iheopiveteni DiBinoloh, of ighiotddsg .covela. pigst. LD p : 
— eingbfoot etl ometosnoModtgse tds. eetiigae Tau} mott snigas) 2ctiieagt {S) 


 Pessibly forming mixed erystah mageq airi-ni bodiz.esob ton arg sesddi can 
| —e.AtiryindBexictsiSeniboe Seettnetan 2 atnuoms aninsycitiw eotinoibonga) (8). 4 
ne esub-ophitic augite partially uralitieac €ssdtiaRi1d aisydqonsig 927830) {h)., stati 
- lash Motoe sgeden lbmiguaend. vitaderg, aliaerD oitydqonety bonisig-sni4 (2) ) 7 
7 bronsite ia ahown by bustihte pecidosecyis »T9qsq ei. ai ‘bedigab _ jon at sh 


racenigeaingenghaomerinaeebey Hub A» sodyb yrudgrog oinydaonstg gatbisioni 
-. - Japge area of magnetite in the ot: AMG ydjazodwaslo bodeildng od of ai 
Thinthprisiaaig, ombeibac pod-tsipit “—~ apes. Bi.o39K3 101 z00g, 275 eotmaggxd. 
 — ghiqgktaaictsiaes ate tue-beigsss eaw, youre, :i!oayear A yolley sit Yo dour 3940 
F; , Snedshach bei biQadds disoned noiensixe vied? sos of cele bas ,eovslone olesd ant 
sitongaR ied, de nehy ip 8 betgligqan, yous Wopas, oii to _nolistorqueial ofT 
ie tieolsinieademeiiT -xoiqrmoD sigittA od! yniariot esqyi-doot edt to asitsaqorg 
esinodeag: Westaseronnigt! Iii a2ao-a} polieyitenyal odi bebaagxe ody, ML od 
hoifegietivad deers word fadi: eradius.Jaio{ aii yc pazsdiene, aiff .eakeog od}. to 
lidalgdiosy: sonsd heamoiieeghm-omnagds) adi. gaiazapaea 4) 26).inciiequi lexgves 
labradorite. The total iren ore content ofb>etwheb od AEA AIN oala_ io. xioigid 
Around some of the magnetite patches, some deuteritic MANROEG, TO saolonet, ae 
the ferro-magnesian mioersi6oiignigases | qo}, betogiae, 19%, engmiivege agye@ 
Pyrit@T beech 2 PE Yee - carsicener OVS DAR oR 
crystels es in No. 765. O81 of aomiosgé........ otingrD sinytigonexs) I - 
—— This spiCbvhmar-t6i - OME MII ON? roid Ber wn <2 BOARDED,» — 
4 aiedud asotsds sent sdt is hr sow noitsnimnxe vilgetgorOg 19% #e0ie nidT 
—penictinstab eiteeqoty oitirystttieds bad.bad. gadno beisesl.orl aatA _.2bem.orew 
— Tagpherptegioatmcessijelts atae aol smodwaitaul oped) me23 sup. orem Rosia x9d)ut 
_ teptaoniof eesrititenod cried? wii: moc 1ort-<bine>: namaste oitgetaatied: tan ti ose at 
oe “vou cskiteagond titacyers-pamath.balssgeus .eonm din, sods oi ylinsivolug 
— beagainsie taldopesnonchar duisen wont, ieomn do amuian AUoMBROTIIOS sodter 93 01 
 Thoationbtypaidinte woe: bobiloig sos asd, aoilagisegy et, telial ods gotia_os oaile. enoxt 
— Ssteqaon: soudtiegn pensgrcomibor deaned & 2h GitT AT. osntw sient aim 
jamgieotigeencotrs) ost j Stimpenbi-erise ga man (ii olodidgme bas stigus jinobsidsi to 
od} didmopsiabiqaseonn do -aoi!gzeace isigr ibovitirglsio. bag bopiviieus ous elsragin 
- eprint ot-qu ediat misrbedse ci-oliscbasdel odT ...besijinueeuse wlisool ots eingalst 
ich eterna gtialego- dane, 0 AeA Rin gt aa) OE -banaiat_ bas bongs anol 
— aftive,,.paisold, bes ekeckinsed 09973,.03 , hagrade 


— etymos  ginoimdq rt ND prs Jjgues.. ot}. om sid, .moitibba al 
= 


The Thermo-Magnetic Properties and History of Some Plutonic Rocks 39 


The Arrigle Complex consists of five major rock-types :— 

(1) Basic Enclaves, gabbroidal to doleritic in texture. 

(2) Tonalites, ranging from true Tonalites through Monzotonalits to Adame- 
lites; these are not described in this paper. 

(3) Granodiorites, with varying amounts of skialiths of sedimentary origin. 

(4) Coarse granophyric Granite. 

(5) Fine-grained granophyric Granite, probably a marginal phase of No. 4; 
it is not described in this paper. 

In addition, there are the usual dyke associates of a plutonic complex, 
including granophyric porphyry dykes. A full description of the Arrigle Complex 
is to be published elsewhere, (by D.J.B.) 

Exposures are poor, for there is a cover of light bog and some glacial drift 
over much of the valley. A magnetic survey was carried out to assist in mapping 
the basic enclaves, and also to trace their extension beneath the Old Red Sandstone. 
The interpretation of the magnetic survey necessitated a study of the magnetic 
properties of the rock-types forming the Arrigle Complex. This was carried out 
by H.M., who expanded the investigation to cover the thermo-magnetic properties 
of the rocks. It is considered by the joint authors that from these investigations 
several important facts concerning the thermo-magnetic, and hence geological, 
history of the rocks may be deducted. 

2. Petrology of Specimens 
Seven specimens were selected for examination : 


4,- Basic. Enclaves... ists. des Specimens Nos. 765, 768, 769, and 770. 
1. Granophyric Granite........Specimen No. 180. 
2. Granodiorites ..............Specimens Nos. 134 and 136. 


Thin slices for petrographic examination were cut at the time the rock cubes 
were made. After the heated cubes had had their magnetic properties determined, 
further slices were cut from the actual cubes. Both sets of slices were examined 
to see if any petrographic changes could be seen in their constituent minerals, 
particularily in those with known or suspected thermo-magnetic properties. Due 
to the rather heterogenous nature of most of the rocks, and consequent changes 
from slice to slice, the latter investigation has not yielded any striking results. 

Basic Enclave No. 765. This is a normal medium-coarse gabbro, composed 
of labradorite, augite and amphibole, with magnetite-ilmenite ; the ferro-magnesian 
minerals are uralitised and chloritised with seperation of rare epidote, while the 
felspars are locally saussuritised. The labradorite is euhedral, in laths up to 10mm 
long, zoned and twinned. The augite is interstitial to sub-ophitic, and is much 
changed to green hornblende and chlorite with granular epidote and fibrous 
actinolite. Some brown hornblende is probably primary. There are some chlorite 
areas composed of chrysotile laths forming a regular pattern but lacking any lime 
mineral. The iron ore forms about 3% of the rock. It is mainly in laths of 
ilmenite up to 1mm long, occuring mainly in the ferro-magnesian minerals. A 
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little fine-grained magnetite is seen in some of the uralite areas, while skeleton 
magnetite occurs in the primary bown hornlende. No leucoxene was seen. Pyrite 
is present in small amounts, simulating magnetite in its mode of occurrence and 
possibly forming mixed crystals. 

Basic Enclave No. 768. This is a coarse gabbro, with basic labradorite and 
sub-ophitic augite partially uralitised. There is some primary brown hornblende 
and some green hornblende which may also be primary. The former presence of 
bronzite is shown by bastite pseudomorphs; it is even possible that the rock once 
contained primary olivine. Rods of haematite may be seen in some of the felspar. 
In general the magnetite occurs in patches, up to 0.5mm long, but there is one 
large area of magnetite in the slide which is 2.5mm across. There is also a 
a sprinkling of fine secondary magnetite formed during the uralitisation and 
chloritisation of the original ferro-magnesian minerals. Pyrite is common, often 
in quite large crystals. 

This specimen was oriented in the field, and so the cubes cut from it have 
aiso been oriented. 

Basic Enclave No. 769. This coarse gabbro closely resembles No. 769, but for 
the fact that there is a great increase in the number of haematite rods in the 
labradorite ; such heamatite is an exsolution phenomenon and is not unusual in 
labradorite. The total iron ore content of the rock lies between 4% and 5%. 
Around some of the magnetite patches, some deuteritic biotite has developed from 
the ferro-magnesian minerals; in part, it must replace some of the magnetite. 
Pyrite is common, but in small grains or agglomerations of grains, not in large 
crystals as in No. 768. 

This specimen was oriented in the field, and so the cubes cut from it have 
also been oriented. 

Basic Enclave No. 770. This is a porphyretic dolerite, with a scattering of 
large plagioclase crystal set in a fine-grained base with a somewhat gabbroidal 
texture. The large felspars are altered to tremolite-actinolite and zoisite: if it 
were originally labradorite, it contained no haematite rods. The augite of the 
base has followed the usual augite-to-hornblende and hornblende-to-chlorite changes. 
The total iron ore present is but 2%, much less than in the other basic enclaves. 
Skeletal octahedra of magnetite are seen, and there is ilmenite in rods up to 0.7mm 
long. There is a very little pyrite present. 

Granophyric Granite, No. 180. This is a representative of the late cross- 
cutting granites which are the most homogeneous rocks of the Arrigle Complex. 
It is a leuco-granite, with quartz, orthoclase and albite forming 97.3% of the 
rock. The albite is subhedral and comparatively fresh. The orthoclase tends to 
be poekilitic; it is fresh but some areas are clouded with a dusting of fine 
haematite, apparently introduced from the Old Red Sandstone which till recently 
covered the outcrop. The quartz is strikingly granophyric; intesgnartie are 
mainly with the orthoclase, but occasionally the albite is involved, suggesting a 
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ternary eutectic (see Nockolds, 1947). Green biotite forms 2% of the rock. It is 


strongly pleochroic, with X a pale buff and Y and Z a clear dark green; it is 


almost uniaxial and the birefringnence is low for a mica. It varies in size from 


well-formed flakes up to 2mm long to clots of numerous small flakes; the centre 
of such clots is usually sieved with quartz granules in a symplectite growth with 
the biotite. This suggests that it is secondary after a silica-richer ferro-magnesian 
such as hornblende or amphibole. The lime released in such a change now 
appears as epidote, grading towards clinozoisite and into allanite. Zoned epidote 
is seen, and pleochroic haloes occur in the biotite around the epidote/allanite. No 
zircon was seen. There is a little reddish-brown and grey sphene : in small grains 
such sphene is almost indistinguishable from allanite. In a few places, alteration 
of the biotite to chlorite has resulted in the formation of sagenitic rutile. There 
is a little pyrite present, pseudomorphed by goethite/lepidocrocite. There is 
a very little magnetite in the biotite, but the total iron ore content of the granite 
does not exceed 0.5%. 

After heating, the granophyric quartz remain unchanged. The pink colour 
of the orthoclase, and the haematite dust, have been emphasised. But the amount 
of ore in the felspar has not increased, nor has the haematite changed to 
magnetite. The albite is clouded; sericite, sometimes stained yellow, has developed 
and also clinozoisite/zoisite, but no chlorite or epidote are seen in the altered 
plagioclase. The yellow staining may be goethite or lepidocrocite, and this would 
involve changes in the iron present. Of all the minerals, the mica shows the 
most marked changes. It is now red-brown in colour, pleochroic from brownish 
yellow to dark foxy red. The epidote appears darker in colour. Magnetite may 
be seen in the mica; it appears to be greater than that noted before heating, but 
the difference is not greater than might be expected from one slice to another. 
The haematite in the felspar has not changed to magnetite, and there is no reason 
why the haematite in the mica should behave difterently. No change was observed 
in the small amounts of allanite and sphene. 

Granodiorite No. 134. This specimen represents a comparatively homogeneous 
portion of the granodiorite, but in it ‘skialiths’ (* shadowstones,’ Goodspeed, 1948, 
p. 67) may still be recognised; these consist of biotite hornfels and patches of 
granular quartz. In the body of the rock, the dominant mineral is oligoclase, 
usually heavily sericitised. The perthitic orthoclase is pockiloblastic and com- 
peratively fresh. Biotite forms about 5% of the rock and shows a very patchy 
distribution. Small poekilolithic biotites may be seen in the skialiths while clots 
of coarser biotite occur throughout the more homogenised body of the rock. 
Where enclosed on orthoclase, the biotite is fresh; elsewhere it is much altered 
to chlorite, in which some granular epidote and sphene occur. Zircon and apatite 
are minor accessories. There is almost no magnetite visible in the slices cut for 
examination. 


Granodiorite No. 136. This is a representative of skialith-rich granodiorite, 
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and in the hand-specimen slivers of biotite hornfels and quartz blebs are visible. 
The slices show the quartz blebs to consist of aggregates of rounded grains 
cemented by a thin coating of limonite. The hornfels consist of small_biotites in 
a base of quartz and clouded oligoclase. There is some poekiloblastic orthoclase 
in the body of the rock, which closely resembles that of No. 134. The biotite is 
much altered to chlorite, in which sagenitic rutile is common ; one small crystal 
of goethite was seen with the rutile. The amount of magnetite is extremely small, 
but there are a few specks of pyrite, with a goethite/limonite rim. The rock is 
so heterogeneous that each slice will vary considerably from the average for 
the rock. 

3. Laboratory Procedure 

In order to make the magnetic observations, the specimens had to be cut 
into cubes with sides 2cm long. This was done by means of a specially designed 
rock-cutting machine designed by H.M. in co-operation with Dr. Bruckshaw. It 
consists of a diamond-impregnated copper disk spun at 3,000 rpm by an electric 
motor. Where oriented cubes were required, the initial cut was made parallel to 
the paint-marks on the field specimen, the operator holding the rock in his hands. 
Thereafter the remaining five faces were cut with the aid of a jig, and a true 
cube cotained from the randomly-shaped specimen ; if the latter had been oriented 
in the field, the resulting cube or cubes were also oriented. In all 25 cubes were 
cut from the 7 field samples examined. 

The method of making the magnetic observations and calculations has been 
fully described by Bruckshaw and Robertson (1948); the same apparatus and 
procedure were used in this investigation. The direction and intensity of the 
remanent magnetisation was determined for each cube in three perpendicular 
directions. From this the amount and direction (for oriented specimens), of the 
remanent magnetisation were determined. The susceptibility was also found for 
each cube by measuring its behavior as ‘an alternating magnet under the action 
of the changing energizing field’. 

These observations were made on the cubes before and after they had been 
heated to above 600°C in a furnace. The results are tabulated in Tables 1 and 2. 

No attempt was made to determine the Curie point directly. By heating 
samples of Nos. 134 and 180 to 300°, 500°, and 600°C successively, it was determined 
that the Curie point for these rocks was higher than 500°C, and probably close to 
600°C: see Chevallier et Piere, 1932. Specimen No. 134 showed an appreciable 
increase in susceptibility (13) and a large increase in remanent magnetisation (32) 
after heating to 600°C ; whatever mineral altered to magnetite, ‘or another magnetic 
mineral (?), did so close to the Curie point of magnetite, which lies around 570°C 
for pure magnetite. Specimen 180 showed almost nil (1.3) increase in susceptibility 
but a marked increase (45) in remanent magnetisation ; whatever mineral already 
in the rock took-on its latent residual magnetisation on cooling around 600°C, did 
so close to the Curie point of magnetite. Similar results were obtained by Michel- 
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Levy et Grenat, (1929), Grenat (1930), Nagata (1943) and others in similar investiga- 
tions into rocks. It is concluded that the magnetic constituent is magnetite. 

Due to the low concentration of the magnetite in the natural rocks and to 
its low susceptibility, no attempt was made to extract the amount required for an 
X-ray powder analysis. 


4. Magnetic Properties 
The following magnetic properties were determined for the rocks in their 


natural state :—The susceptibility (&,) and the remanent magnetisation (#) for all 
specimens and also the angle and direction of magnetic dip for the oriented 
specimens. With its remanent magnetisation, a rock possesses an induced magnet- 
isation due to the earth’s magnetic field (H); the value of this induced magnetisa- 
tion, k,,H, is readily obtained by multiplying the susceptibility by the intensity of 
the earth’s local magnetic field. In the laboratory, H equals 0.545 oersted. The 
ratio of R to k,H is called ‘7, and is a measure of the magnitude of the 
remanent/induced magnetisation. 

Table 1 shows the results obtained. At least two cubes had been cut from 
each rock specimen, and in all 25 cubes were investigated. Where approximately 
the same figures were obtained on two or more cubes cut from the same rock, 
the values are not repeated in Table 1, but suffixes ‘a,’ ‘b,’ etc. have been added 


to indicate the number of cubes measured. 
Table 1. Magnetic properties of natural rocks prior to heating. The general accuracy equals 


10 10-® c.g.s. units. Susceptibility and Remanent Magnetisation are expressed in 
10-* c.g.s. units. N.D. means ‘not determind’. 


‘ . wi gae teak ak es Pe 
Rock Type and Susceptibility | Remanent Mag.| Ratio R/k,H . 
Specimen No. | kn R 7 | r . ae on ye | 
Basic Enclaves wn ~ oe vonelernase i) 

765, a, b 105 120 2.0 ;| =- | = 

768, a | ND. | 370 N.D. | —33° E. 148° | 

768, b N.D. 260 N.D. —25° E. 162° ! 

768, ¢c 270 270 | 1.8 - cee 

769, a 375 810 | 4.0 —38° E. 115° 

769, b 155 525 6.2 —30° EB. 115° 

770, a, b, c 90 40 | 0.8 = a / 
Granites | 

180;°a,"B; ce, “a «| 125 <10 <0.15 _ — | 
Granodiorites ; 

134, a, b, c, d | 20 <10 <0.9 - -— 

136,: a, .b, c,d, | Nil Nil - wa 

e, f, g - 
—* i _— = | 


5. Thermo-Magnetic Properties 


When the foregoing measurements had been completed, the cubes were heated 
to 600°C in a furnace and then allowed to cool slowly to room temperature in the 
earth’s magnetic field. In so cooling through its Curie Point (which could not be 
greater than 580°C), each cube acquired new magnetic properties. 
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These new magnetic properties were measured in exactly the same manner as 
that used on the unheated-treated rocks. Table 2 shows the results obtained. 
The new susceptibility is indicated by ‘k’ and the new remanent magnetisation 
by ‘P’; this may be described as the thermo-remanent magnetisation of the 
specimen. The ratio of P to &,H is indicated by ‘Q’. This ratio is of considerable 
interest and diagnostic use. It was first introduced by Koenigsberger (1933, 1938) ;- 
Thellier (1938) and Nagata (1941, 1943) have shown that it is a magnetic character- 
istic of rocks, independant of the field strength for small field intensities. 


Table 2. Magnetic properties of the rocks after they had been heated to 
600°C, and allowed to cool in the earth’s magnetic field. 
Susceptibility and Remanent Magnetism are expressed in 10- 


c.g.s. units. 
Rock Type and | Susceptibility | Remanent Mag.| Ratio P/kxH | yp, 
Specimen No. ke a Q »P 
Basic Enclaves | 
765, a 165 320 3.5 | GE? 
765, b 165 330 2.8 59° 
768, a 350 1,225 6.5 = 
768, b 270 540 4.3 — 
768, c 290 1,520 9.7 64° 
769, b 285 1,150 tae 64° 
770, a 390 1,025 4.8 68° 
770, b 310 680 4.0 65° 
770, ¢ 345 840 4.5 69° 
Granites | 
180, a 230 420 3.4 65° | 
180, b 135 235 3.2 68° | 
180, c 155 290 3.5 67° 
Granodiorites | 
134, a 330 570 3.2 68° | 
134, b 115 200 ase 67° 
134, c 335 570 3.8 63° 
136, a 400 830 3.8 67° 
136, b 375 735 3.6 64° 
136, c 275 620 4.1 68° 
136, d 250 470 3s5 67° 
136, e 275 600 4.0 TO” 
136, f 330 685 3.8 69° ; 


It will be noticed that the heated cubes take on the magnetic dip approxi- 
mating to that measured for the laboratory, (65°); the mean value for the cubes 
is 66°, and the range extends from 59° to 70°. Sample 769, a completely disinteg- 
rated during heating, and 769, b was severely fractured. This phenomenon is 
rare and may be due to the expulsion of volatiles or to the release of strain at 


moderately high temperature. 
The ratios of the three sets of values given in Tables 1 and 2 have been 
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calculated and are given in Table 3. The ratio of the susceptibility after heating 
to the susceptibility bofore heating is called ‘X’; any departure of this ratio from 
unity may be taken as implying a formation or destruction of magnetic constituents 
in the rock by and during the heat treatment process. The ratio of the thermo- 
remanent magnetism after heating to the remanent magnetism of the natural 
rock is called ‘S’; its interpretation is by no means simple. Finally, the ratio of 
Q to T is called ‘Z’; this Z ratio is the S ratio corrected for changes of suscepti- 


bility produced by the heat treatment. 


Table 3. Ratios between the Susceptibilities, Remanent-Induced Magneti- 
sation and Remanent Magnetisation after and before heating 


to 600°C. 
Rock Type and Ratio k:/Rn Ratio P/R Ratio Q/T ) 
Specimen No. Xx po Ss 4 ) 
i: Basic Enclaves | 

765, a 1.6 2.7 ) 1.7 | 

765, b 1.6 2.8 . 1.8 

768, a N.D. pa Bee. N.D. 

768, b N.D. ~ By N.D. 

768, ¢ isi 3.6 5.4 

769, b 1.8 2.0 1.2 

770, a 4.3 26. 6.0 

770, b 3.5 17. 5.0 

WO c 3.8 a1. 5.6 
Granites ) 

180, a i rg 40. 23. 

180, b ee 25. 16. 

180, c ya 30. a 
Granodiorites 

134, a 16. 55. 3.5 

134, b 6. 20. 3.5 

134, ¢c 1%. 55. 3:5 

136) a 40. large - ; 

136, b 40. large _ 

136, c 30. large 2 

Mao, 25, large - 

136, e 30% large des 

186,7 £ 35. large én 


6. Discussion 6.-1. Magnetic Properties 

The figures for remanent magnetisation given in Table 1 indicate a sharp 
division of the rocks into two groups; a similar, but less sharp, division may be 
made on the figures for susceptibility. The basic enclaves show appreciable 
susceptibility and remanent magnetisation; the granite and granodiorites are 
markedbly low susceptibility and almost no remanent magnetisation. This general 
grouping of the rocks investigated corresponds to the known geological facts: 
basic rocks are normally more magnetic than those at the acid end of the series. 
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Considering the basic enclaves as a group, it is seen that the three gabbroes 
(Nos. 765, 768 and 769) are very similar. Susceptibilities range from 375 to 105, 
remanent magnetisation from 810 to 120, and T values run from 6.2 to 1:8. . The 
magnetic properties of the fourth enclave, No. 770, are very different, and the 
rock is a porphyretic dolerite. Susceptibility is 90, the remanent magnetisation is 
but 40, while T is 0.8. Other workers in this field have tentatively established 
the criterion that rocks with T values above 3, and possibly above 2, have 
crystallised and cooled through their Curie point, (Nagata, 1943). This criterion 
would establish the high-temperature origin of the gabbroes and at the same time 
suggest that the porphyretic dolerite was of low-temperature origin. 

The granophyric granite shows appreciable susceptibility, and a very little 
magnetite was seen in the slices. Hence the extremely low remanent magnetisation 
indicates that the rock has not cooled from above its Curie point, or that new 
low-temperature magnetic minerals have been added after cooling, or that the 
rock has an abnormally low coercivity. The last explaination is rejected, for 
magnetic rocks in general have a high coercivity. The thermo-magnetic properties 
will permit the elimination of one of the two remaining explanations. 

The granodiorites have low to nil susceptibility and remanent magnetisation. 
Virtually no magnetite was seen in the slices. The lack of susceptibility also 
indicates that no magnetic minerals are now present; hence little can be learnt of 
the thermal or geological history of the rock from an investigation confined to its 
magnetic properties. 

Information regarding the absolute direction and sign of the dip of the 
basic enclaves is based on two oriented specimens obtained from in situ outcrops ; 
from these oriented cubes were cut. The error lies within 5°. Table 1 shows 
that the dip is around 30°, and that there was a reversal of magnetism. The 
declination for one specimen was 155°, and for the other 115°. These results do 
not necessarily indicate a reversal of the earth’s magnetic field towards the close 
of the Caledonian orogeny. The basic enclaves were presumably free to move and 
rotate within the fluid granodiorite, whether the latter is considered as a magmatic 
or metasomatic product. Clearly this linc of investigation would enable the 
amount of relative movement undergone by different enclaves to be determined, 
but only significantly if a great number of in situ enclaves could thus be investi- 
gated. 

6. Discussion 6.-2. Thermo-magnetic Properties 

Basic Enclaves. The 765 samples have increased their susceptibility by 60% 
and their remanent magnetisation by a factor of 2.7, giving a Z ratio of 1.7. 
This Z value implies that the rock has lost 40% of its remanent magnetism if the 
earth’s magnetic field was then of the same strength as it is today, or that the 
minimum field in which the rock cooled was 0.31 gauss. The 769 b sample 
increased susceptibility by 80% and doubled its remanent magnetisation, giving a 
Z ratio of 1.2. This Z value implies that the rock has lost 17% of its remanent 
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magnetisation if the earth’s magnetic field was then of the same strength as it is 
today, or that the minimum field in which the rock cooled was 0.46 gauss. The 
differences in the estimates of the minimum terrestrial field some 400 million years 
ago may be due to different coercivities in the two rocks, or to the fact that only 
minimum values can be estimated. 

Only incomplete results are available for sample 768. The high Z ratio 
would suggest that the rock has not passed through its Curie point, but this is 
not supported by other evidence. 

The 770 specimens have increased susceptibility by 29022, and their remanent 
magnetisation by a factor of 21, giving a Z ratio of 5.5. The increase in suscepti- 
bility is sufficiently large to indicate that some mineral which can thermally change 
to magnetite has been deposited in the rock after cooling, or that the rock has 
never before been heated above its Curie point. The high Z ratio also implies 
that this rock has never before been above its Curie point, for there are no 
reasons to believe that it has demagnetised more rapidly than the other rocks of 
the group. 

Granite. The 180 samples have increased their susceptibility by the low 
factor of 30%, but their remanent magnetisation by the extremely large factor of 
32, giving a Z ratio of 21. Even allowing for the fact that the initial remanent 
magnetisation may have been in error by 50%, since it was so low, this change is 
far greater than that found for other rocks in this present investigation, or in the 
literature on this subject. Likewise the Z ratio is the highest encountered, but 
the Q ratio is very close to the average for all the rocks from the Arrigle 
Complex. 

Granodiorites. The 134 samples have increased their susceptibility by a mean 
factor of 13 and their remanent magnetisation by 36; the resultant Q and Z ratios 
are quite normal at 3.2 and 3.5 respectively. The 136 samples have behaved in a 
similar manner, though the increased susceptibility shows a factor of 33 and the 
remanent magnetisation factor is of the order of 50; however the Q ratio is quite 
normal at 3.8, and the indeterminate S ratio may also be normal. The heat 
treatment has therefore brought about the formation of new magnetic minerals 
in the rock, as witnessed by the high X values. As the rock cooled, these new 
magnetic minerals passed through their Curie point, and the great increase in 
the thermo-remanent magnetisation is shown by the high S values. The new 
magnetic mineral is taken as magnetite. 

6. Discussion 6.-3. Thermal and Geological History of the Rocks 

The magnetic changes produced by heating a rock above its Curie point and 
allowing it to cool in a magnetic field whose strength is close to that of the 
earth’s may be considered under two heads :— 


(i) The heating may alter the magnetic properties of the magnetic mineral 
already present in the rock. 


(ii) The heating may cause new magnetic minerals to form from the pre- 
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existing minerals presesent in the rock. 

The effect of (i) should be to alter the remanent magnetisation yet leave 
the susceptibility unchanged or very little changed. The effect of (ii) should cause 
the susceptibility to increase, and since new magnetic minerals will have been 
formed, the remanent magnetisation will also increase. Of course the two eftects 
may be produced in the same rock undertreatment, while other rocks will contain 
no magnetic minerals before or after the heat-treatment. Such peculiar effects as 
the destruction of magnetic minerals by the heat-treatment, or the imparting of 
‘latent magnetic properties’ to pre-existing minerals are not considered. And it 
is assumed that the rocks will have the normal high coercivity. 


Table 4. Utalisation of the natural magnetic and the thermo-magnetic ratios of rocks 
to assist in the diagnosis of their thermal history. 
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It may be accepted that magnetite is the magnetic mineral with which this 
investigation is dealing. It has been indicated (Herroun and Hallimond, 1943, 
p. 220) that the only satisfactory explanation for the remanent magnetisation 
found in many natural rocks is that the magnetite took on its permanent remanent 
magnetisation in cooling through its Curie point. The Curie point for pure 
magnetite has been determined as lying between 580° and 570°C ; impurities lower 
this temperature, (Chevallier et Piere, 1932). The Curie point for the rocks from 
the Arrigle Complex has been fixed roughly as around 550° to 580°C (see Section 
3 above). 

Some temperatures at which silicate minerals crystailise are of interest for 
comparison. The albite-oligoclase plagioclase felspars commence to crystallise 
around 1400°C. The eutectic for quartz-orthoclase is 800°C and the hypertectic is 
680°C. Quartz-orthoclase-hornblende has a ternary eutectic point at 680°C. Primary 
magnetite in such rocks is usually considered to be one of the first minerals to 
crystallise ; it should do so at temperatures high above its Curie point, and under 
such circumstances should take-up a strong remanent magnetisation. When it 
does not do so, the theory that it has cooled from a high-temperature ‘magma’ 
falls under suspicion. 

The application of these arguements to the basic enclaves show that the 
three gabbroes (Nos. 765, 768 and 769) are normal igneous rock which cooled from 
above the Curie point of their magnetic constituents. They fall into Case A of 
Table 4. This is in harmony with the petrography and presumed geological 
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history of these rocks. 

The porphyretic dolerite (No. 770) occurs as a basic enclave in the grano- 
diorite, but appears to have a very different thermal history from the gabbroes ; 
it falls into Case D of Table 4. Magnetite is developed on heating to above the 
Curie point, and all the indications are that the rock as at present constituted 
has not previously passed through its Curie point. There are two possible geologi- 
cal explanations of the phenomenon in this case. After the dolerite had cooled 
through its Curie point, low-temperature metamorphic changes may have destroyed 
all the magnetic minerals; these reformed when the rock was again heated in this 
investigation ; the metamorphic changes may have occurred when the dolerite was 
engulphed by the granodiorite. It is not easy to see why the same changes did 
not affect the gabbroes. On the other hand, the dolerite might be a metamorphic 
rock, derived say from volcanic ash, passing through a hornfels stage into a 
pseudo-dolerite. The porphyries rather favour this theory, but otherwise the 
petrography denies it. The investigation has focussed attention on the uncertainty 
of the origin of what would otherwise be taken for a normal igneous dolerites. 

The granophyric granite (No. 180) falls into Case B of Table 4; the suscepti- 
bility has barely increased while the remanent magnetisation has increased enor- 
mously. This indicates that the rock, as at present constituted, has not previously 
passed through its Curie point. At first sight, this appears as a startling result, 
for the geological evidence points to the granophyric granite as a late cross-cutting 
intrusive granite, while the petrography indicated that it has crystallised from a 
molten silicate smelt. The iron ore introduced from the Old Red Sandstone (or 
even from some other source) does not affect the question. It is seen to be 
haematite before heating and to remain haematite after heating. However, the 
petrographic examination of the slices indicated that there had been late-stage 
deuteritic alteration of hornblende into biotite with exsolution of silica as quartz 
blebs. At the same time, the magnetite seen in small quantities in the biotite 
may have formed ; such magnetite would be low-temperature magnetite and would 
have low or nil remanent magnetisation. On heating in the laboratory, the existing 
magnetite would pass through its Curie point and then obtain its remanent (or 
thermo-remanent) magnetisation. This is considered to be the correct explanation. 
The main body of the rock crystallised at high temperature; the magnetite is 
deuteritic, but its idioblastic habit makes it liable to be identified as primary 
idiomorphic magnetite. 

Both the granodiorites (Nos. 134 and 136) fall into Case D of Table 4; thus 
they resemble the porphyretic dolerite No. 770). New magnetite appears to have 
been created; it could form from low-temperature minerals belonging to the 
skialithic portion of the granodiorite, or from low-temperature or decomposed 
high-temperature minerals belonging to the body of the granodiorite. From this 
it is clear that the magnetic, or potentially magnetite forming, minerals of the 
granodiorite as found in the field, are low-temperature minerals, whether or not 
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they originated from high-temperature minerals. The examination did not reveal 

what minerals were changed into magnetite on heating; Michel-Levy (1929) con- 

sidered that sphene decomposed on heating, but it could hardly form appreciable 
amounts of magnetite. > 

It is concluded that the relation between the magnetic and the thermo- — 
magnetic properties of rocks yields much information as to the thermal and 
geological history of the rocks, but that this information must be used with caution 
and must be governed by the history of the rock as revealed by its petrography 
and other geological features. 
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The Diurnal Variation of Cosmic Rays 


By Kazuo NAGASHIMA 
Earth Science, Kyoto University 


Abstract 


An attempt is made to examine the energy gain or loss of 
the primary cosmic ray produced by an electro-magnetic field of 
solar stream, and then to explain the diurnal variation of cosmic 
rays at the time of the magnetic storm by this energy variation 
of cosmic rays. 
The amplitudes and the times of the maximum intensity of 
the diurnal variations are calculated at various altitudes and lati- 
tudes by using the above mentioned electric field. 
Whether the diurnal variation of cosmic rays averaged over 
all days can be explained or not from the standpoint of the 
electric field hypothesis is also discussed. 
1. Energy variation of Cosmic Rays in an Electro-magnetic Field Produced by a 
Solar Stream 
It has been supposed that the occurences of the magnetic storm and aurora 
were due to the solar stream which was composed of neutrally ionized particles, 
emitted from the active region of the sun, with a mean velocity of about 
2-10%cm/sec., [1 and 2]. Although the outbreak mechanism of the stream from the 
sun is still obscure, it is sufficient to quote the experiments of Meinel [3] and 
Gartlein [4] in order to show its existence. Such a stream will be polarized 
electrically by the solar magnetic dipole field, and moreover if the electrical 
conductivity of the stream is very good, the magnetic field in the vicinity of the 
sun will be frozen in the stream and pushed forward with it [5]. Although the 
influence of the electro-magnetic field of the stream upon the diurnal variation of 
cosmic ray intensity was pointed out by several authors [6, 7, 8, 9 and 10], its 
quantitative discussion is not yet given. As the electrical influence of the stream, 
polarized electrically by the solar magnetic dipole field, upon the diurnal variation 
of cosmic rays is of the same character as that of the stream, freezing the 
magnetic field in itself, so hereafter only the latter case is considered. Strictly 
speaking, the magnetic field frozen in the stream is not uniform but, in order to 
simplify the analysis and to obtain the essential feature of the influence of the 
stream upon the diurnal variation, the uniform magnetic field is assumed to be 
frozen in. First, we will consider the following case. 
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A stream of Equatorial type—This denotes the one which is bounded by the 
two parallel infinite planes, perpendicular to the solar equatorial plane. The 
magnetic field H frozen in it is parallel to the boundary planes and moves with 
a velocity v. On the other hand, it is assumed for simplicity that the magnetic 
field outside the stream is zero. Such a situation is shown in Fig. 1. As seen 


yy SUN from the coordinate system (x, y, z) 


fixed to the earth (hereafter refered 
to as the rest system), an electric 
field is observed besides the magnetic 
one, and then the cosmic rays change 
their energies by this electric field 
when they encounter the stream. In 
order to examine the cosmic ray orbit 
and its energy variation in such an 
electro-magnetic field, it is useful to 

Fig. 1. A stream of Equatorial type. introduce another coordinate system 
(x’, »’, 2’) fixed to the stream (hereafter refered to as the moving system), because 
there exists only the magnetic field in this system. When a cosmic ray, having 
its initial velocity uw, and its initial total energy W, as seen from the rest system, 
reaches the boundary surface on the left hand side (cf. Fig. 1), we make the 
rest system transform into the moving one by the Lorentz transformation as 


follows, 
rn — (Uo). s/1—F? +) — (Uo)yt+v + — (lo), /1— #2 
Gy ata EE, wt) oe (11) 
145Gb "14 pltie 1+ptey 
where p= : é 


and a prime indicates the quantity measured in the moving system. In this system, 
the cosmic ray is subjected to only the magnetic deflection and reaches another 
boundary surface with a velocity u’. The relation between wu’ and u', is given by 


the following equation, 


’ , a 
cos¢ =COS¢? o— 
p” (1.2) 


As for the notations in Eq. (1.2) confer Fig. 1. Here again we make the 
Lorentz transformation and obtain the final velocity u and the final total energy 


W of the cosmic ray in the rest system. 


=~ u', J/1-# u,= swt ’ u,= —— ° (1.3) 
ft 3 / : u a 
1p» — ces 


The relations of the final quantities to the initial ones in the rest system are 
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obtained from Eqs. (1.1), (1.2) and (1.3), and the results are as follows, 


4W, = W—W,=feaH, (1.4) 

uy 4W, 
( € eos Ao ates dee B W, (1 5) 

cosd cos¢y= P 57 - 
| (to “49 4W, 4 G3 

J ( c ) W, W, 
. JWi(mey 7 
sin A= JV Wo+ 4 OMe sin ho, ( ) 

in which Uy = Uy 


and a is the width of the stream. 4W, in Eq. (1.4) is the energy gain of cosmic 
ray, after traversing the stream from the left to the right, and is the same for 
all particles regardless of their energies and velocity directions, if they have the 
ability of passing through. On the other hand, the particles, passing through 
from the right to the left, lose their energies by the amount of 4W,. The 
minimum energy W,”(4) of the particle, which has the ability of passing through 
the stream and the final velocity of which is directed towards 4, is approximately 
given by the following formula, 


Ww,” 
W""(A)= a (1.7) 
7 Aa 
where 1 age ee yi (1.8) 
The above relations hold under the next conditions. 
3 4W, 
B<l1, Bacosa> (1.9) 


W,” denotes the minimum energy in case that A=0. If we adopt W,” as an energy 
unit and represent the total energy by the next equation, 


W=XW,", (1.10) 


a dimensionless variable X can be adopted as a new measure of the energy. 

As seen from Eqs. (1.5) and (1.6), the particles after traversing the stream 
are directed towards the given directions (4, ¢), which form a definite region in 
the velocity space. Next, we examine this region under the conditions of Eq. (1.9) 
which result also in the following relations by considering Eqs. (1.6) and (1.7), 


4W,<W (or W,) } 


sinA~sin A). (1.11) 


Case [I] when the observer is located at a point P, (x=a). 


(1) Xcosi=1 


In such a case, at least some of the particles, incident upon the stream from 
the left with the given energies, can traverse it as seen from Egs. (1.7), (1.8) and 
(1.10). By considering Eqs. (1.5), (1.6) and (1.11), their velocity directions at P, in 
the rest system are restricted by the following relation, 


K. NaGasuiMa 


Fe T=P>(Ya) os (Accelerated Region) (1.12) 
Cae cos(1 -' eae - 4,18) 

This holds approximately under the conditions of Eq. (1.9) or (1.11) provided that 
a cosa> f. (1.14) 


Such a situation is illustrated in Fig. 2a. The 
particles, their velocity directions lying in the 
striped region of a circle in Fig. 2a, are those, 
(Ge which have passed through the stream and then 
* Ac, Reg. gained their energies by the amount of 4W,. On 
the other hand, those, having their velocities lying 


in the remaining region of the circle, do not gain 


or lose their energies. This situation is shown in 


Fig. 2b. It is the purpose to examine the diurnal 
AWe 


GO Wr rw 2m(p) 
Fig. 2 a. (upper) and b. (lower) ‘ . : a . 
Velocity distribution, and 4W, in Fourier series regarding 4W, as a function 
energy variation 4W as a func- of ¢, (cf. Fig. 2b). 
tion of 9, in case that x=a and 
X cos?.=1. 


variation of cosmic rays produced by such an energy 
variation mechanism, so that it is useful to expand 


AWAg)= +a cos¢+0,sin¢)+(a,cos2¢+b,sin2¢)+---- 

= 9 + Aicos(y—6,)+ Ascos(2—6z)+ +--+. (1.15) 
a in the above equation contributes to the intensity variation not dependent on 
local time. This contribution is somewhat different from the one in the case of 
the static electric field [9] because a is in general a function of W and A, but is 
not refered to here. The diurnal variation of cosmic rays is caused by the 
quantities a, and 5, or A, and 4, which are obtained easily from Eqs. (1.12) and 
(1.13) as follows, 


: 24W, ; A 
a= — AM sino.).= a ee (1-49 ') ‘ 


(1.16) 
24W, seca 
py ate [1+cos(¢.)-J= t (:- Xx ), 
AW, o> aetaly sec 
A, = * /Al+coGalel = gf 1-~yX (1. 16)’ 


tand,= — a/ 4 cosA—1 e 
Hereafter —a,, 5; and A, are refered to as the 12 hour, 18 hour and resultant 
components of the diurnal variation of 4W, respectively, because the vectors —a, 
b, and A; point to the earth from the directions of 12, 18 and 7, (=2x—d,) hour 
local times respectively. It must be emphasized here that although all the 
particles traversing the stream gain (or lose) their energies by the same amount 
of 4W,, A: depends on X and 4 as seen from Eq. (1.16y. 
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(2) Xcosa<1. 

In this case none of the particles can traverse the stream so that we are 
not able to observe the accelerated particles at P,; and then a,, 6, and A; are zero. 
Case [II] when the observer is located at a point P in the stream (0<x<a). 

(1) Xcoss=1. 

Similarly as in the above case, the particles passing through the left boundary 

" surface gain their energies at P by the amount of 4W{=fexH) and are restricted 
in the velocity directions as follows, 


(¢..=¢=(¢.)-+ (Accelerated Region) (1.17) 


On the other hand, particles, traversing from the right, lose their energies at P 
by the amount of 4W,-_,. (=fe(a—x)H) and are restricted in the following velocity 
directions, 

27 4+-(¥,)-—=9—[¢-]es (Decelerated Region) (1.18) 


where 


x or) 


(¢.).= cos'( 1-5 78 (z>(v.).=0) (1.19) 


[e.].= cos-*(—14+ 5% .?804), (on >[0,].=n) (1.20) 


Such a situation is illustrated in Fig. 3. a, and 5, or A; and 4d; in Eq. (1.15) are 
obtained as follows, 


a= AW et sin{y.1.—sin(¢.).] a PAM | 


us 


b, = A Wat costy,), —cos{¢,].]= =e i *(1- — 


Ai = AW J 2(1-sin{¢,]-sin(¢,), —cos[¢,].cos(¢,).], 
cos(¢?,).— cos[¢,],. 
Bans sin[¢,].—sin(¢,), * 


(1.21Y 


In particular, when the observer is located at the centre of the stream (x=a/2), 
these quantities become as follows, 


a-—24Wa [seci secd 
= 


~ ¥ 2—-—y—}> 
_ 24W, secA 1.22 
b= (i 7) (1.22) 
2 
Ay= AW, => +(es).. 


ria ’ 


It is a remarkable fact that as seen from Eq. (1.21), the 18 hour component (5;) is 


independent on x and only the 12 hour component (—a,) varies with the movement 
of the observer’s position. 


(2) Xcosé<1. 
In the same way as above, the velocity distribution is divided into the 
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following three regions, 
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Velocity distribution, and energy 
variation 4W as a function of 9, in 
case that 0<x<a@ and Xcos)=1. 


Xcos)=1. 
22—(G2)-—=9—(Pe)e» (Accelerated Region) 
[¢v.].>¢>2z2—(¢,)-, (Forbidden Region) < 
An—[¢.].—¢veI[¢.]., (Decelerated Region) he 
and 2n+(¢.)->¢>4n—[¢,].. (Forbidden Region) 


where (¢,). and [¢,]. are defined by Eqs. (1.19) and (1.20) respectively, (cf. Fig. 4). 
The Accelerated and Decelerated Regions have the same meanings as in the 
previous case. The Forbidden Region signifies that particles having their velocity 
directions lying in this region do not exist. This region is caused by the action 
of the magnetic field which is called the magnetic cut-off. So that in this energy 
region where Xcosi<1, the contribution of the magnetic field to the diurnal 
variation of cosmic rays coexists with that of the electric field. It must be 
emphasized here that both contributions have only the 12 or 24 hour component 
and not the 18 or 6 hour component, because the distribution of the above three 
regions as a function of ¢ is symmetrical about the point g=z (or 2z) as seen from 
Eq. (1.23) or Fig. 4. The contribution of the electric field can be obtained, in the 
same way as in the previous cases, by assuming that 4W is zero in the Forbidden 


Region. 
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a= = {—4W,sin(¢,).+4Wa-.sinly,1-3 | (1.24) 
bd, =() 


On the other hand, the cotribution of the magnetic field is obtaincd as follows: 
Let i(X)dX be the cosmic ray intensity in the energy interval between X and X+dX. 
If the influence of the electric field upon #(X) is neglected, i(X) is constant in all 
the regions except for the Forbidden Region in which 7(X)=0, and is expanded in 
Fourier series regarding as a function of gy. The 12 hour component —a,' of «X) 


is given by the following equation. 
a,'=— seein {sin(¢,).+sin[¢.]-}, (1.25) 
6,'=6. 


The term a,‘/2 produced by the action of the magnetic cut-off contributes to the 
cosmic ray intensity variation not dependent on local time, but is not refered to 
here. 

As above, we have examined the influence of the solar stream of Equatorial 
type upon the diurnal variation of cosmic ray. It is noteworthy here that, when 
the observer goes into the stream from the outside, only the 12 hour component 
—a, (or —a,‘) is changed in its magnitude, on the contrary the 18 hour component 
b, (or b,*) is always constant irrespective of the observer’s position. In case that 
Xcosi=1, this remarkable character is 
illustrated by the harmonic dial in Fig. 5 
in which the radii of large and small 
circles O and O' are 24W,/x and 
4W,/nx respectively. If we choose P on 
the Ob,-axis so as to satisfy the following 
relation OP=24W,(1—secd/X)/x and 
draw a straight line parallel to Oa,-axis 
through P then PP; and OP represent 
—a, and b, respectively, when the ob- 
server is located outside the stream, (cf. 
Eq. (1.16)). The observer entering in the 
stream, the 12 hour component —a, 
begins to increase and continues to do 12"(-a,) 
so until the observer reaches the centre Fig. 5. 
of the stream, where —a, is represented 
by PP;, as seen from Eq. (1.18) or (1.19). In case that Xcosi<1, such a character 
is easily seen from Eggs. (1.24) and (1.25). It is also noteworthy that even if the 
magnetic field vector is directed inversely, the above character is unchanged except 
for the point that the 18 hour component is taken the place of by the 6 hour 
component. 


In the above discussion, the uniformity of the magnetic field is assumed in 
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the stream. If there would be some deviation from the uniform magnetic field, a, 
and b, given in the previous equations would be changed of course. But. even. so, 
this deviation from the uniform magetic field would not change essentially the 
above mentioned character that only the 12 hour component is changed in its 
magnitude when the observer goes into the stream from the outside. The defor- 
mation of the stream from the Equatorial type also changes a, and 3,, but would 
not alter decisively the above character. At any rate, it will not be meaningless 
to examine the influence of the deformation upon a, and }. 

Fig. 6b shows a stream with a finite cross section in xz-plane and with 
infinite length along oy-axis. The magnetic field is uniform inside the stream and 
its direction is parallel to oz-axis. At the upper and lower boundaries of the 
stream in Fig. 6b, the magnetic lines of force change their directions and are 
parallel to oy-axis. As this magnetic field, its direction being parallel to oy-axis, 
does not produce any energy variation of cosmic ray because the direction of this 
magnetic field is parallel to that of the stream’s velocity, so the influence of this 
field is neglected here. By the same procedure as discussed above, a calculation 
of the energy variation of cosmic rays is made in a case when the observer is 
located at the centre of the stream which has the relation b=a. Only the result 
is shown in Fig. 6 which shows some deviation of A; and 7, with the increase of 
the laiitude from those of the stream of Equatorial type given by Eq. (1.22). As 
discussed above the magnetic field at the upper and lower boundaries does not 
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Fig. 6 a. A; as a function of X, in case that b=a. Numbers attached to curves indicate 
latitudes. 
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Fig. 6 b. 7; as a function of X in case that b=a, where T;=2x—&,. Numbers attached 
to curves indicate latitudes. 


cause any energy variation of cosmic rays, but produces the deflection of cosmic 

ray orbit, so that the neglection of this magnetic field may produce some error in 

T, at the high latitude. If 5 becomes very longer than 4g, it is sufficient to consider 

only the influence of the stream of Equatorial type. 

2. Diurnal Variation of Cosmic Rays produced by the Electric Field of Soiar 
Stream 

In the preceding section, we have obtained the 12 hour and 18 hour com- 
ponents (or the resultant component) of the energy variation of cosmic ray 
produced by the electro-magnetic field of solar stream. Now we examine the 
diurnal variation of cosmic ray in the earth’s atmosphere expected from the 
existence of such 12 and 18 hour components. 

As pointed out in a previous paper [11], the intensity variation of cosmic 
ray caused by an electric field not derived from a potential is represented as 
follows: In the normal state when any disturbing field except for the earth’s 
magnetic field does not exist, the vertical intensity of cosmic rays N(A, p) observed 
under the atmospheric depth pg/cm® at the geomagnetic latitude A degree, is given 


by the following equation, under the assumption that all the primary cosmic rays 
consist of protons. 


NG, x)=\"m(E, EAE. 2.1) 


Where m(E, p) and «(E) are the over-all multiplicity [11 and 12] and the differential 
intensity spectrum of the primary cosmic rays respectively, and E, represents the 
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cut-off energy at 4°. When an electric field not derived from a potential is added 
to the earth’s magnetic field, the intensity varies as follows, 


Ni, p, 4E(E, A, t), SE} = \ mG pu +L{E-4E(E, i, 1), 4E(E, 4, t)¥] 


xi{E—AE(E, 4, )}dE (2.2) 


4E(E, 24, t) is an energy gain which the primary cosmic ray, having its kinetic 
energy E at 4° on the earth, has undergone during its passage through the electric 
field, and is in general a function of E, A and local time ¢ because JE depends 
upon the cosmic rays trajectory. L{E, 4E} is Liouville’s effect [11] and given by 
the following formula. 


Loe an. Ae ee 


(E+ mioc?)?—(myc?)? * (2.3) 


6E, in Eq. (2.2) is the variation of cut-off energy produced by the electric field, 
but is neglected hereafter under the assumption that its influence upon the intensity 
variation is small [13]. Then the intensity variation relative to the normal state 
is defined by the following formula. 

AN _ N{A, p, 4E(E, 4, )}—N(A, D) 

N N(A, p) 
In the case of the above discussed electric field of solar stream, the formula of 
4E(E, t) is determined under the following assumption that only the earth’s 
magnetic field is considered in the region where its influence upon the motion of 
cosmic ray is superior to that of the solar stream, on the contrary, outside this 
region only the field of solar stream is taken into consideration. The deflection of 
cosmic ray orbit in the earth’s magnetic field is obtained from the results of the 
model experiments of Malmfors [14] and Brunberg [15] and shown schematically 
in Fig. 7, where 4y represents the latitude of the asymptotic orbit of cosmic ray 


(2.4) 


Fig. 7. 


and ¢, “is {the difference between theYlongitude of the asymptotic orbit and that 
of the observed point. BN in Fig. 7 represents the direction from which a cosmic 
ray is incident into the earth’s magnetic field with a maximum energy gain, and 
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is hereafter refered to as the direction of anisotropy. If the equatorial plane of 
the earth is assumed for simplicity to be paralled to that of the sun, this direction 
coincides with that of the vector —a,, b, or Ai. Then 4E(E, 4,4) is given by the 
following formula, 
AE(E, A, t)=4E\(E, dy) cos (Yrt+t—T), (2.5) 
where 

—a,(E, Ay)=4Wa-fa(E, Ay) 
AEE, hy)=4Warf(Es Ay)= 4 OE, dy) = 4 Wear fr Ey dy) i 
AE, A4y)=4WarfalE, Ay), ; 

T\(E, Ay)=2n — 6,(E, Ay). 


If the resultant component A, is adopted as 4E,, T, in Eq. (2.5) is a function of 
E and dy as seen from Eqs. (1.16)’, (1.21) and (1.22) or Fig. 6. On the other hand, 
if the 12 hour component —a, (or the 18 hour component b,) is adopted as 4E,, 
T, is independent on E and dy, and is equal to 2z (or 3z/2). Introducing Eq. (2.5) 
into Eq. (2.4), we are able to obtain the diurnal variation of cosmic rays by the 
following method: the values of 4N/N are calculated numerically at the twelve 
local times, ¢=0", 2", 4",....22%, respectively, and then the diurnal amplitude (4N/N), 
and the time of the maximum intensity fu.x are obtained respectively by the 
harmonic analysis of the above twelve values of 4N/N. As an example, we will 
calculate (4N/N), and fix by adopting A; and 7, in Fig. 6. In the following calcula- 
tion, Neher spectrum [12 and 16] given by Eq. (2.7) is adopted as #(E) and m(E, p) 
obtained in the previous paper [11] is used. 


1B) perc SO DEVI (E=0.8 Bev.) Te 
=0. (E<0.8 Bev.) 
If it is assumed that 
4W,=1.3-10'ev. and = 5-10", (2.8) 
it follows from Eq. (1.5) that 
W,"=10" ev., (2.9) 


which corresponds to X=1. In such a case, (4N/N);, and fmax of the ionizing and 
neutron components at various altitudes and latitudes are shown in Table L If 


we determine an energy E by using the value of fmx in Table I so as to satisfy 
the following relation, 


Prttmx—T,=0, (2.10) 


we see that the cosmic rays, whose kinetic energies are FE, being incident into the 
earth’s magnetic field from the direction of anisotropy, just hit the point of the 
geomagnetic latitude 4 degree on the earth at the local time imax. We are able to 
define this energy E as the mean energy of the diurnal variation of cosmic rays. 
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ionizing neutron _ 
ae : _|___ component _compnent _ 
Ps/em") 1930 | 0 680 - 
i AN \ Jaen 
Ce) | 0.30 | 0.78 | 0.59 
Table I. i ee 100 16 6 56m!) gh 33m 
Diurnal amplitude (AN/N);, time of the maxi- EBe. | 9 | 19 | be 
| 
mam intensity fmax, mean energy E and S in ies f 5 
aang of the resultant component A; shown in - Me Boe ical 
ig. 6. ‘AN 
: _— Gc. 0.64 | 2.0 17 
(AW, =1.3-108 ev. and}.8=—— -10-*) oe 
tea 11» 14m 8h 20™ 8b 42m 
50°| — pitas 
E Bev. 21 7.0 8.2 
5 0.26 0.30 0.29 
~) FANY 
hoe), %| 0.36 | 2.8 
tmax 136'36™ | 125 32m 
80° —— 
E Bev. 13 9.0 
S 0.09 0.52 
Table II. 


Diurnal amplitude (AN/N),, time of the maximum intensity tmax, mean energy E and S. 


(A) denotes the case of the resultant component A; shown in Fig. 6. 
(a) Se the 12 hour component —a, 
(b) ~ the 18 hour component 5, 


(AW. =0.6-10% ev. and $= 5 -10-) 
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ionizing neutron | inoizing neutron ionizing neutron 
a component | component . component | component | component | component 


~ee/e™) 1030 | = 680 | 1030 | 680 1030 | 680 
(| (az % 0.1 | 0.32 | 0.10 0.22 0.16 0.31 
| tmax | 11816" gh 32m 541m 5b15m | 13b37m | 12h19m 
; EBev. | 25 22 ‘ z 
s | 0.14 0.27 : it bal scat Sim Sn OP" 
Ca % | 0.35 1.0 | 0.25 | 0.99 0.30 0.51 
r 12> 31 gh 56m Bh 4gm sak ee 756m | 1526m | 14n54m 
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E, thus defined, is shown in Table I, together with (4N/N), and tnx. 


In case where 


4W,,=0.6-10' ev., a=% 10-2 and W,"=5-10" ev., (2.11) 


(4N/N):, tmax and E are given in Table II. In it, the diurnal variations caused by 
the 12 hour component —a, and the 18 hour component 4, in Fig. 6 are also 
tabulated respectively. 

In order to find the relation between (4N/N); and E, it is useful to simplify 
Eq. (2.4) as follows, 


an _ Sg QE AEE An) C08 (Cott 19) mE, PEE 


Me [", m(E, p)i(E)dE 


=4q(E): 4E,(E, Ay) 608 pti 7) (2.12) 
under the assumption that 7(E)=(E+mpc*)*. gE) in Eq. (2.12) is given by the next 
equation. 
“ 2 
Qé)= 1-11 +E/mey~° (2.13) 
A bar in the above equation denotes the mean value with respect to the energy. 
If it is supposed that the diurnal variation is produced by the primary cosmic ray, 
whose energy is equal to the mean energy E obtained in the above, its amplitude 
becomes 


ie - 4W, 


E+ moc 

as seen from Eq. (2.12). (4N/N); in Table I or II does not necessarily coincide 

with this amplitude but can be connected with this by the following relation, 

( 
N 


QE Ey = flE, dy)-qE)- 


4 - =, 4W. 
) a9 AE iy) Eg eS aE) ae ary, (2.14) 
where 4W, denotes the energy gain of the cosmic rays which traverses the stream 
and s is called the smoothing factor [21] which is caused by the deflection of oe 
ray in the earth’s magnetic field and f is called the effiiciency factor which is 
caused by the influence of the stream’s form upon 4E(E, Ay). If it is assumed 
that 4E does not depend upon E and dy and is equal to 4 W., it follows that the 
cosmic ray being incident upon the earth loses its energy by the amount of 4W, 
regardless of its velocity direction and then that the field which produces such i 


character of 4E must be a static electric field. In 
such a case, Eq. (2. 
written as follows [17], Eq. (2.14) can be 


4N = ke) AW, | 
NS eR me (2.15) 


From the above consideration, S may be interpreted as the deviation of the cosmic 
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ray variation, produced by an electric field not derived from potential, from the 
variation produced by a static electric field. By using the value of (4N/N), and E 
in Table I and II, S is obtained and tabulated in the same place. As seen from 
the comparison between Tables I and II, the values of E and S depend upon 4W, 
to a certain degree, even if the same type of the stream is adopted, but as their 
dependences on 4W, is very weak they can be regarded as nearly constant. 

3. Discussion and Conclusion 

1) The electric field of solar stream. 

As seen from Tables I and II, it is necessary for the value of 4W, to be 
equal to 10*ev. in order to produce the diurnal variation of the vertical intensity 
of the ionizing component at sea level, its amplitude being about 0.2 or 0.5 percent. 
If it is assumed that the velocity of the stream is 2-10° cm/sec. and the width of 
the stream is 5-10'cm, the magnetic field H, which causes such an energy variation 
(4W,~10° ev.), is obtained from Eg. (1.1) as follows, 


H= aes =10-° gauB. 
The solar magnetic dipole field, its strength at the sun’s surface being about 1 
gauf at the equator, becomes about 1.3-10-’ gaug in the vicinity of the earth’s orbit 
and is much smaller than the above required value. Then we are not able to 
expect the diurnal variation from the solar magnetic dipole field. But, as already 
discussed in paragraph 1, if it is supposed that the magnetic field in the vicinity of the 
active region of the sun is frozen in the stream and pushed froward with it because 
of the good electrical conductivity of the stream [5], the difficulty can be avoided. 
2) The diurnal variation of cosmic rays at the time of the cosmic ray storm 
(Ds). 

It is a well-known phenomenon that the intensity of cosmic rays is decreased 
at the time of the magnetic storm (Dst) and the anomalous diurnal variation is 
observed at the same time (Ds), [18]. This phenomenon is called the cosmic ray 
storm [19]. Ds is defined as the difference of the diurnal variation vector of the 
disturbing day from that of the quiet day. Various authors [6, 7, 19 and 20] 
examined the Ds; phenomenon and pointed out its many characteristics. rere 
[21] analyzed Ds by using data at many stations and concluded that the direction 
of anisotropy which causes the observed Ds is on an average towards about 13 hour 
local time. As pointed out in paragraph 1, when the observer goes —_ the aegis 
from the outside, only the 12 hour component —a, is increased in its sananities> 
and the 18 hour component 6, remains constant, (cf. Fig. 5). 1 the cosmic ray 
storm is assumed to begin at the time when the earth just enters into the stream, 
then the direction of the newly added anisotrapy vector is towards 12 hour local 


time and its magnitude is 
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at the centre of the stream as seen from Eqs. (1.16) and (1.22), (cf. Fig. 5). This 
newly added 12 hour component given by Eq. (3.1) is shown in Fig. 8, as a function 
of Xcosd. Such a theoretically expected direction of the newly added anisotropy 
vector (12hour local time) seems to coincide with the observed one (13 hour local time) 
above mentioned. It must be emphasized here that as pointed out in paragraph 1, 
even if the magnetic field vector in the stream is directed inversely, the direction 
of anistropy does not change at all. Yoshida compared the energy dependence of 
the newly added 12 hour component shown in Fig. 8 with the one analytically 
obtained from various observations. 
This comparison seems to show that 
both dependences coincide with 
mutually in the order of magnitude 
[21]. 

Although the main character 95 
of Ds as above seems to be explained 


=, 


by this mechanism, the fine structures 
of Ds, for example, the phase ad- 
vancement of the anisotropy of the 
Ds vector in the course of a cosmic 0 


(7, 19 d 21] 1 10 XCOS A 10 
ray, stormys Ctes.. sees , an 

. . ‘ Fig. 8. The amplitude of the newly added 12 hour 
cannot be explained by such a simple component as a function of Xcoa), at 
model of the stream. But, these the centre of the stream. 


problems may be explained by a suitable choice of the stream’s structure. 

3) The diurnal variation of cosmic rays averaged over all days. 

Although there is enough ground for controversy as to whether this averaged 
diurnal variation is nothing but Ds or contains some components caused by other 
mechanisms quite different from that of Ds [22], what is examined here is whether 
the observed variation can be approximately explained or not by the following 
assumption: this averaged diurnal variation is caused by an electric field not 
derived from a potential, (hereafter refered to as the electric field hypothesis). 
The diurnal variations calculated in paragraph 2 are those caused by the electric field 
of the stream, hence it may be difficult to compare them directly with the observed 
averaged diurnal variations. But, as it seems that even other form of the 
electric field, which is different from that of the stream, does not lead to any 
result quite different from those obtained in Tables I and II as long as the electric 
field hypothesis is accepted, so it may be possible to discuss qualitatively this 
problem according to the results in Tables I and II. Another difficulty of making 
a direct comparison between the observed and theoretical diurnal variations is as 
follows: the theoretically calculated diurnal variations are those of vertical inten- 
sity, on the other hand the observed variations are caused by the cosmic rays, 
incident upon the instruments with finite apertures. Taking these difficulties into 
consideration, we examine the observed averaged diurnal variation from a stand- 
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point of the electric field hypothesis. 

The energy variation 4W, of 10°ev. can produce the diurnal variation of 
about 0.2 or 0.5 percent at sea level as seen from the above tables, on the contrary 
in the case of the static electric field, the same amount of 4W, can produce the 
intensity variation of about 1 percent or more at sea level [11 and 17). Then it 
may be concluded that the diurnal variation is less sensitive to 4 W, than the 
intensity variation. This less sensitivity is caused by the factors f and s discussed 
in paragraph 2. The mean energy E in the same table agrees nearly with the 
observed one [23], although the former depends upon the stream’s form. 

Altitude dependence.—As seen from Table I, the ratio of (4N/N), of the 
ionizing component at the upper atmosphere to that of sea level is equal to 2.6 at 
the geomagnetic equator and 3.1 at the geomagnetic latitude 50°. If the amplitude 
of the averaged diurnal variation at sea level is assumed to be equal to about 0.2 or 
0.3 percent which is a reasonable value, the expected amplitude at the upper atmos- 
phere (p=0) is between 0.5 and 1 percent. This does not contradict the observed fact 
that the amplitude of diurnal variation is less than about 1 percent at the altitude of 
balloon flight [24]. On the contrary, as discussed in the previous paper [11], the 
magnetic field hypothesis of the diurnal variation [25] or others, in which primary 
cosmic rays in the comparatively low energy region play an important role in the 
diurn2! variation, will produce a larger altitude dependence than that of the electric 
field hypothesis because low energy rays are apt to be absorbed in the earth’s 
atmosphere. The underground experiment of Mac-Anuff [26] shows that at the 
time of the large magnetic storm when the intensity decrease became about 5 
percent on the ground, the intensity of underground (60 m.w.e.) did not change, on 
the contrary the underground amplitude of diurnal variation increased as first 
pointed out by Yoshida [21]. The theoretical altitude dependence of (4N/N), is 
smaller than that of the intensity variation in the case of the static electric field 
as seen from Table I, and Fig.7 in the previous paper [11]. So that this character 
of these altitude dependences can explain the above cited underground experiment, 
although this is a problem of Ds and must be discussed in 2). The time of the 
maximum intensity fmax also changes with altitude. At a constant latitude, fmx at 
the upper atmosphere is in general two or three hours earlier than that of sea 
level as shown in Tables I and II, but it happens that the former is later than 
the latter in special cases. 

Component dependence--The observed ratio of (4N/N), of the neutron 
component at Climax (A4~50°N) to that of the ionizing component at Freiburg 
(A=50°N) is between 3.2 and 3.8 [15]. As pointed out by Yoshida and Kamiya [27], 
the intensity variation of the ionizing component at Freiburg is very small 
compared with the ones at other stations which are located at nearly the same 
latitude as Freiburg. It seems to be more adequate to adopt the intensity variations 
at the latter stations, then the values of the observed ratio cited above become 
small. The theoretical ratio is equal to about 2.6 in the case of Table I or 2.9 in 
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the case of Table IJ, and then seems to be in agreement with the observed one, 
although the theoretical one depends more or less on the form of the electric 
field. 

Latitude dependence—The theoretical ratio of (4N/N), of the ionizing (or 
neutron) component at geomagnetic latitude 50° to that of the ionizing (or neutron) 
component at geomagnetic equator is about twice the observed ratio of the 
corresponding component [28 and 29]. Such a discrepancy between the theoretical 
and observed values may be due to the following facts: firstly the theoretical one 
is that of the vertical intensity, on the contrary, the observed one is that of the 
omnidirectional intensity; secondly the theoretical one depends more or less upon 
the form of 4E(E,A,?) given by Eq. (2.5). fax of the ionizing component at low 
latitude (A=0°) is one or three hours earlier than the one at high latitude (4 = 50°). 
This tendency is caused principally by the deflection of cosmic ray in the earth’s 
magnetic field and seems to be in agreement with observation [29 and 30]. 

As above, the electric field hypothesis seems to be suitable to explain also 
the averaged diurnal variation except for the latitude dependence, although it is 
somewhat difficult to discuss quantitatively this variation on account of the above 
mentioned causes. 

In conclusion, the diurnal variation of cosmic rays at the time of the cosmic 
ray storm may be explained by the energy variation mechanism of cosmic rays 
in the electro-magnetic field of solar stream. The fine structures of the observed 
Ds cannot be explained by such a simple model of the stream as discussed in 
paragraph 1, but may be explained by a suitable choice of the stream’s ftructure. 

The necessary energy variation 4W, is about 10° ev. in order to produce the 
diurnal variation of the ionizing component at sea level, its amplitude being about 
0.2 or 0.5 percent. This requires that the magnetic field is about 10-° gau8 in 
the vicinity of the earth’s orbit. Such a requirement cannot be satisfied’ by the 
solar magnetic dipole field, but may be achieved by assuming that the magnetic 
field in the vicinity of the active region of the sun is frozen in and pushed forward 
with the stream. 

Some parts of the averaged diurnal variation of cosmic rays also seems to 
be explained by the electric field hypothesis. 
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Meeting of the Society of Terrestrial Magnetism and Electricity : 
The 16th General Meeting was held at the Nagoya University on November 3-5, 
1954. 
Number of the Reports read at the Meeting: 
‘Cosmic Rays, 7; Atmospheric Electricity and Atmospherics, 8 ; Rock Magne- 
tism and Electricity, 8; Atmospherics and Radio Meteorology, 9; Geomagne- 
tism, 11; Ionosphere, 10. 
Tanakadate Prize was awarded for the following excellent workers: 
The 16th, Mr. K. Sinno, Mr. T. Obayashi and Mr. H. Kamiyama ; 
Study on Ionospheric Storms with Special Reference to Disturbance-Daily- 
Variation. 


The 17th General Meeting was held at the Tokyo University on May 6-8, 1955. 
Number of the Reports read at the Meeting; 
Atmospheric Electricity, Radio Meteorology and Atmospherics, 11; Cosmic 
Rays, 13; Rock Magnetism and Electricity, 13; Geomagnetism and Earth 
Current, 11; Ionosphere, 14. 
Tanakadate Prize was awarded for the following excellent worker: 
The 17th, Mr. N. Kawai; 
Studies on Stability of Remanent Magnetisation of Sedimentary Rocks. 
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